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INTRODUCTORY REMARKS 



In his- memoir Helraholtz (reference 1) showed the . 
possilDility of matheciat ical analysis of those types of 
flov/ of incompressihle liq^uids that are 'characterized 
by the formation of so— called rays (Strahle) or jets 
within the region at rest. Following the work of 
Helmholtz a rather large nuin*ber of investigations, 
devoted to the same problem^ appeared in foreign and 
Russian scientific literature. At present the fully 
worked out Joukowsky method (reference 2). permits the 
solution of any problem on steady, irrotational flow 
ofan ideal liquid under the following conditions: 
first, the fluid throughout moves parallel to a certain 
plane^ the- flow being bounded by plane walls perpendic- 
ular to this plane, and secondly, the motion takes place 
in the absence of external forces. (The same conditions 
are imposed in almost all problems of this type,) 

The analogous problem for an ideal gas has hardly 
been touched upon. The author is familiar with only one 
paper which deals with gas j^ts; namely, the one by P. 
Molenbroek (reference 3)/ Kolenbroek set up the differ- 
ential equations, on which the problem of gas jet flow 
depends and gave certain particular integrals of these 
equations; these equations, however, hardly correspond 
even to the theoretically conceived .motion of the gas. 

In the present paper a method is presented with the 
aid of vrhich it is possible, in many cases, to find the 
solution of a given problem on the fjow of an ideal gas. 



* Scientific . Memoirs, koscow University, 1903, 
pp, 1^121,. 
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Here it is necessary to impose all the conditions which 
are assumed in the analogous problems on an incompressible 
liquid as previously mentioned by Joukowsky (reference 2), 
But, in addition, the applicability of the analysis here 
developed is further restricted by the special req.uirement 
that the velocity of the gas particles must nowhere ex-- 
ceed the velocity of sound for the particular physical 
state of the gas at a given point (local velocity of 
sound)^ Corresponding restrictions are likewise imposed 
on the limits within which the pressure may vary. If 
this additional condition is not satisfied, stable 'motion 
apparently is not possible. It is assumed, however, that 
with the aid of a certain hypothesis, stated in this 
paper, the problem can be analyzed also for the case 
where the additional condition is not satisfied. The 
mathematical treatment of this problem, however, is left 
to another paper, 

A brief summary of the contents of this paper is 
presented here^ 

In part I the differential equations of the problem 
of a gas flovr in two dimensions is derived and the partic— 
ular integrals by which the problem on jets is solved are 
given. Use is made of the same independent variables as 
Molenbroek used, but it is found to be more suitable to 
consider other functions, The stream function and veloc- 
ity potential corresponding to the problem are given in 
the form. of series. 

The investigation of the convergence of these series 
in connection with certain properties of the functions 
entering them forms the subject of part II. 

In part III the problem of the outflow of a gas from 
an infinite vessel with plane walls is solved. 

In part IV the impact of a gas jet on a plate is 
considered and the limiting case v/here the jet expands 
to infinity changing into a gas flow is taken up in more 
detail, Jhis also solved the equivalent problem of the 
resistance of a gaseous- medium to the motion of a plate. 

Finally, in part V, an approximate method is pre^ 
sented that permits a simpler solution of the problem of 
jet flovrs in the case where the velocities of the gas 
(velocities of the particles in the gas) are not very 
large. 
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A number of supplementary notes are appended at the 
end of this report, the second of v/hich establishes a 
relation between the analysis of part V v;ith certain 
problems, in the theory of minimal surfaces, 

: A further interesting remark may be noted here: 
The\results obtained in parts III and IV, at least 
quali tat iyely^ agree sufficiently vrell v/ith test results, 
although the experimental investigation of the phenomena 
accompanying the jet f crma t i on was conducted under con- 
ditions very different from those assumed in these theo- 
retical inves t igat i ons . 

The principles of the method with its application 
to flows presented here were briefly communicated to 
the Moscow kathematical Society at the beginning of 1896, 
A more detailed presentation was made at the eleventh 
session of the experimental scientists and doctors in 
1901. 

In concluding these introductory remarks deep 
appreciation is expressed to E. A. Bolotov for his kind 
help in proofreading the manuscript. 
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PART I 



GENERAL PRINCIPLES 01 THE METHOD OP INVESTIGATION 



An infinite mass of a perfect gas contained between 
two parallel planes is assumed . and, in addition, "bounded 
"by certain cylindrical surfaces perpendicular to these 
planes. One of the latter is assumed as the coordinate 
plane XY^ Let the gas be in stabilized motion and let 
the direction of the velocity throughout be parallel to 
XY, The effect of external forces will be neglected and 
it will be assumed that the velocities have a potential^ 
Since it is desirable to avoid vorticity formation, it 
is necessary to consider the pressure as a function of 
the density^ It is convenient to take 



and thus assume an adiabatic process. 

The magnitude Y equal, for atmospheric air, to 
1.4025 (reference 4) is the ratio of the specific heats. 
It is preferred to consider the motion as constant heat . 
process in view of the small heat conductivity and radia- 
tion of the gas particles. Because of this the adiabatic 
process at large velocities appears most closely approach- 
ing the true conditions. In any case, the result of this 
analysis must be considered as a first approximation for 
the reason that no account is taken of the connecting 
chains between the jjarticles and the resulting viscosity 
forces, friction at the walls, and so forth, factors 
which, in the case of -gas flows, are possibly of greater 
effect than in the case of liquid flows. 

Under the foregoing assumptions the velocity potential 
cp is a function of x and y and, for the components of 
the velocity u, v, the expressions 



p « kp 



(1) 



u 



Sep 

11' 




by 



(2) 



With the density of the gas denoted by p the con- 
dition of continuity is written 
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^e^. \E1^ 0 • (3) 
OX by ' 



The Bernoulli law in this problem may, with the aid 
of equation (l), he reduced to the relation 



;2\ P'Y 

P = Ro ' ' 



where 



= u^ + v^, a - —yZT^* ^ " " constant ( 



the density at the point of the gas is evidently p 
where V = 0, ^ 

Per hriefness, 

l! = T 
2a 



so that 

Equation (3) indicates the existence of a function 
\l/ determined hy the equations 



po Po 



Prom equations (2). and (6) with the aid of equation 
(5) a relation is obtained between the functions cp and 
\1/ given by the formulas 
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(7) 



= - (1-T) 



ocp 

1y 



The function ^i^ represents the stream function, 
tho equation ^ = constant "being the equation of a stream- 
line. By assigning successive constant values and 

Cg to the latter, it is readily shown that (C^- Cg)po 
expresses the mass of gas per second flowing through a 
cross section of the jet "between the streamlines. 



Equations (?) are transformed by taking cp , for 
the independent variables and considering x, y as 
functions of cp and The relations 



Sep by* &cp bx 

b\|/ by* b^i^ bx 

are readily obtained where 

bcp b\!; bcp b\l/ . v ^ s ' 
bx by by bx . 

From the foregoing equat ions the reciprocal of the square 
of the velocity is obtained: 




Equations (?) become 
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■ (8) 

It Is necessary to consider the derivatlTes with 
respect to tp and ^ of the irariablee t « Y^/Sa and 

. 0=arctg|=arctg|=-arctg| (8'); 

the inclination of the velocity to the X = axis is evi- 
dently e. Differentiation of with respect to ^ 
results in 

or, on the hasis of equation (8) 
Differentiation of 8 with respect to <p gires 



or 



4 



These relations lead to the equation 



Further 



do 

whence, with the aid of equation (8), there is obtained 
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But 



4aT» 

By making the final reduction in the formula for and 

substituting t for V® in relation (9) gives the 
relat i ohs 

= 2t ( 1-T ) — 

0\l7 t<p 

Q (10) 

Sep 1-(2P+1)T b^l/ 



Passing next to the independent variables 9 and 

taking cp and ^ as functions of the former gives the 
formulas required 



(11) 



E equat io ns (,11 ) and Xl2^ cons t i tute a solut ion of the 
pr ohlem of the flow of a gas if ihe"*ran^e of variables 
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.•lli^iiB lli^ £l^5£ 2> 1 £2^5:5 "^A^H i^££ik:££ with 

The 9 region will te 'considered singly connected"^ 

and . c-lo se,d. 

In order to show that the function \1/ is fully 
defined for the given conditions it will "be proved that 
the contrary is net true^ Let it he assumed that there 
exist two functions and xl/^ satisfying all these 

conditions. It will 'be shown that ^i— 4^2 = 0, The 
function = ^i^^s everywhere in the given region of 

values T, G is finite and continuous, satisfies ectua-- 
tion (l2), and at the "boundary of the region assumes the 
value zero, Multiply the left side of equation (12) hy 
^dTd9 and integrate vrithin the limits of the t, 9 re- 
gion. If the result of the integration is denoted hjr 1^ 
suhstituting for' ^ there is rtjadily obtained: 

r r / /^^^A^ i^(2g+i)T. 
r ^^'3 i-(2p+i)T 

where the double integral extends over the entire 9 
region, once over its contour. Since on the contour 
"i/^- - 0 the equation I ^= 0 can he true only if the 
douhle integral "becomes zero. Under the ahove—indicated 
conditions, however, the function under the integral sign 
may he .either p-ositive or zero. It is clear that the zero 
value must be taken, and this leads to the equations 

= 0, -^^-3 =r 0 and = constant 0 

as v/as required to he proved; 
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Among the req.uired conditions for the existence of 
a definite solution it has been mentioned that throughout 
the region of gas flow the inequality 



1 

^ ^ 2P+ 1 

must be satisfied. The significance of this requirement 
will be explained. Turning to formulas (4) and (4') for 

T = I- 5 -J— = 0.17 (13) 



gives 



2a 2p+l 



= 23+1 ' 1 + Y 



where 

1+Y • 2e+l' ^ Ho ^ P 
whence 

Y-i 



or, making use of relation (l) gives 

.^'=f'^ (13') 



Thus the restriction imposed on t is equivalent 
to the requirement t^ha^t the yeloc ity of the jga^_nowhere 
exceed t^he velocity; of EXo5a^iti£n'"j2f"*5£jjnd for the 
particular physical conditions at the poinT under con**- 
sideration^ It is supposed that such velocities, at 
least for established flows^ cannot even exist, (See also 
reference 5, and the authors cited by him,) 

The limiting value t « l/(2p+l> establishes also, 
the limits within which the pressure may vary in the region 
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occupied- "by the moving gas mass. Thus, if the variable 
T is everywhere. less than its limiting value, then 




But p7 = Y/(Y-1) = 1 + P; hence 

if it is assumed that V equals approximately 1,40, 

The author turns to the derivation of other very 

important theorems with regard to the motion under con- 
sideration to show, in the first place, that iji^ S^l^^ljt^ 
P-2j-§iJil5j; ^ J 2^11:31^^1^^ S: i!jiJJ£lJ;_2ii ^ f ^h e coordinate s , 

• wj- thin the f JL ow i; e^j. o n , ha jlili^J ^ maximum 
21 ^ minimum. To prove this, it might "be possible to con- 
sider onl7 the following condition. If a point existed at 
which cp had a maximum, there would then have to exist 
ahout it a closed curve on which cp had a constant value 
less than the maximum. In such, case the gas would flow 
through this curve from outside to inside the area bouucled 
"by it. The mass, hounded hy the curve would increase v/ith 
time and the motion could not he steady, By similar con- 
sideration the as sump t i on of a minimum of cp is like wise 
shown to he impossible. But since jbhe theorem on the func- 

J^he fact t_hajt ijt j._s true aljo for the coordinate^ J, 
Z^^^-I^^^ f unct io ns of the independent varj.able^ ^cg 
^nd ^, another proof applicable to^'all these""? unc t ions 
al"so will be given, 

?rom formulas (7), 

bx cx oy by 

At a certain point A of the flow region let cp have 
a maximum (or minimum). About this point take a closed 
curve (0) along which cp maintains a constant value k, less 
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than the maximum (or greater than the minimum). kultiply 
the above equation in cp by F(<p)dxdy and integrate 
over the region hounded by the curve (C). Integration 
by parts yields 



yg B« Dt 
^(k) J j (1-T)$ ^ dy+F(k) / j (1-T 



)3 ^ dx 
6y 



^ [ f THcp) (l-T)^ ^fB^ ^ ixiy ^ 0 (a) 

where the symbol / indicates that the function under 

the integral is the difference of the values ( 1- t") 
at the points B» and B (fig. 1) and similarly for 

r« ^ 

/ . Since, in the case of a maximuiii, cp increases 

D 

in i^assing in\\'ard from the contour (C) at points B, 
B» the result is dx > 0 in passing within the 

region of integration. But since, at the first of these 
points dx > 0 and at the second < 0 for the motion 
along EB ' , then 



and similarly 

D I 

D ■ . . . 

The function r(cp) is chosen so that F(cp) and F'(cp) 
are, everywhere within (C), greater than zero. Turning 
now to the atove-der ived re lat i on (a) it can te seen that 
all its terms are less than zero and therefore impossible, 
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hence also the assnmption of a maxiinum cp, (If the author 
assumed a irinimum, the signs in the substitutions would 
have teen reversed, and the function J chosen so that 
everywhere 'P(cr) and P»(cp) < 0, and again would have 
arrived at the impossibility of (a).* 

The same' cons.iderat ion proves the correctness of the 
derived theorem also, for the other above— men t icned con-^ 
ditidns, formulas (a) being required for functions 
x(cp, ^^ ) and y(cp,\l/ ), * 

¥ith the aid of equation (lO) it. is not difficult 
to prove a similar theorem also for the function and 
therefore th_e velocity of the flow H'^ewi^e cannojb have 
^ ffi«:2ElBHS 1^1 £t i^i^: iables ^cg, a niinj-mum _ 

iJi£ minimum value of T j.^ order to 

prove this the following; equation is constructea on the 
basis of formulas (lO) 

h l-(2^+l)T -3^1 br b (1~-t) It , ^ 
(I-t) — + = 0 (b) 



Assume that there exists in the cp , ^ plane a point 
where T has a maximum or minimum, ?ake, about this 
point, a curve (C) v/ith constant value of r.;ultiply 
the oquaticn for T by a certain function f(T) and 
integrate the left part over the area "bounded bj' the curve 
(C), Integrating by parts yields the relption: 




^ (l-'O („] +1 ^ l-^ J Jf t(T)dc?d\lr=0 




In quite the same manner, as in the above— cons idered 
cases, the impossibility of this relation will be proved. 
It is readily seen, however, that the proof will be valid 

only for the condition: T — within the region of 

" 2p+ 1 ■ 
flow., and this condition has already, been assumed and its 
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physical meaning explained. The case T = 0 is itself 
excluded fromthe range of applicability of the above 
considerations and for the following" reason : At the 
point T s: 0, if this point lies within the flow mass, the 
streamlines meet. It i's readily seen that in this case 
the coordinates x, y cannot he s ingle— valued functions 
of cp, the latter region will he represented, at least, 

by a two-rsheet Riemann surface not assumed in setting up 
the double integrals that figure in these considerations. 
It is easy to show, however, without any formulas, that 
the value t = 0 is the minimum 't. For this, it is 
sufficient to remember that ^ = (u^ + v^)/2a, and, 
since this function is everywhere positive, the value 
zero is actually its minimum. In what follows, only 
such gas flows for which the critical point 1=0 lies 
on the bounding contour of the 9 region and the 

surface of the cp.\l/ region, a single sheet surface will 
be considered. 

By setting up formulas (10) the differential equation: 

b , be b t(i--t) be 

T ( 1-- t) + --^ =0 

Sep Sep bV 1-^(2^ + 1) T 

and applying the above— descr ibed device it is found that 
its 9.1?£jL]i!iX cannot have e i ther a 5ijaxi.mum or a 

ml.nimum. In the same vay the absence -of turning values 
also for the functions cp , ^, of t and 0, if the 
latter are taken as the independent variables is estab- 
lished, For this purpose formulas (ll) must be used, 

From the foregoing theorems proved it is clear that 
in the cp, region there cannot exist closed curves along 
which the functions x, y, T, 9 maintain constant values; 
all such curves must end at the boundary of the region, 
Similiar considerations hold for the 9 region and 

the curves cp constant and \l/ = constant. 

In application only such problems as correspond to 
a 9 region bounded by. concentric circles a.nd the 

straight line segments passing through their centers will 
be kept in mind. The magnidutes T, 0 will be taken as 
the polar co ordinat es . of ^t he points of their region and 
the common center of the boundary curves will be the pole 
of the coordinates. 
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Tot these conditions from the theorem on the im— 
possihility of a maximum or minimum of ^(cp,^) G(cp,^l^) 
it may he concluded that inner points of the ' g 
re§ion correspond to inner points of the cp, ^ region^ 
It raay he noted further that by making use of the al>- 
sencG of a maximum or a minimum of the function y\f (t^ e)^ 
there -can again he obtained the theorem already proved 
on the uniqueness of the fuuction if it is continuous 
within the T 8 region and is given on its houndaries, 
^he series of its. houndar^ Z.§iMl S^Z, IR jg£il^ r a 1 , _ a 1 s o 
J^J^ ^ J;^£^n t i,n UjO us, ^ 

A problem on the flow of a gas will now he considered, 
AssiTme as knov/n the corresporiding contour of the region 
of the variables -T^ 9 satisfying the condition t <' 
l/(2p + l) ; finally y\f on the contour is knov/n. If it 
is i)0ssible, from certain considerations, to conclude 
that the given problem has a solution and if a continuous 
function ^ satisfying the given conditions is found, 
then this function i^rill actually represent the stream 
function, since no othsr is passible.- because of the 
theorem on the uniqueness of the solution of the differ- 
ential equations of the same type as the equations for 
the function (See also reference 5*) 

Side considerations, as to the existence of a solution, 
are not, however, always a priori possible, and such being 
the case, having obtained a function \l/ and through it cp , 
it may be questioned as to whether these particular func-- 
tions give a possible solution of the problem^ In order 
to remove such doubt it is necessary to show each time 
that the formulas for cp and ^ determine T and 9 as 
single-value functions of x and y. 

In order to clarify this point,,, the reasoning will 
be, as follows. ' Let a single-valued* function 9) be 

defined; then from the formula 

5x '^^ bcp ^- bx ■ 
5 T bcp b T. d\i/ 0 T 

and similarly on the basis of relation (8) and formulas 

bx ^ co's 0 • Sy _ sin 9 
. ; bcp . V 2aT Sep -.ySaT 
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there Is obtained 



V2aT — = rr-- cos 9 - ~ sm e (1-t) 
and similar ones, 

. Thus the derivatives of x( "J", 6) and y(T,e) are 
determined n,s s ingle— valued functions of t and 9, 

"sin gle-valued f unct i on a of x,£. But 



and from equations (8) and (ll) the relations 



VP I 



frcLi which it is clear that, if everywhere ^ < l/(2B+l), 
the equation (x,y)/(Tje) 0 is possihle only if both 
partial derivatives of the function V are simultaneously 
zero. This can happen at a singular point of one of the 
curves ^ - constant if such singular point on the curve 
exists. In general, 'to deny the existence of such points 
is impossible, but it can be stated- that curves ^( G) « 
constant will in no- case form a loop, since. closed curves 
\1/ ~ constant ^^rould then exist v/ithin the loop. Hence, the 
branches of our curve, after forcing the singular point, 
v/ill support themselves against the boundary of the region 
somewhat as shown in figure 2» If, however, it is known, 
at least from the conditions of the problem, that all the 
curves ^ « constant issue from the same point of the 
contour Q and a^ain meet at another point of the 

boundary, then the abo ve**ment ioned disposition of the 
curve is impossible and therefore the vanishing of the 
(3!:,y)/(T , e) is likewise impossible. The same is also 
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true in the case v/here, starting from the same point of 
the boundary of the region, the curves ^K'^, 6) = constant 
then divide into pencils, each of them again converging at 
one point. 

An entirely different picture will result if a steady 
gas flow under conditions so that exceeds the limit- 

ing value l/(2p+l) is sought. The Jacohian (cp , ^i/)/ 
( 9) in the region of Q where T is greater than 

the limiting value will then he the difference "between - 
two positive q^uantities and will "become zero along a cer- 
tain curve. Consider, for example, the case where to 
the ho.undaries of the cp,^ region there corresponds in 
the 9 region the semicircle ACB and its diameter AB, 

the center of the semicircle heing at = 0; let ^ 
along this contour have some constant value'^ 

Along OA evidently ( b^/ e) O; on the semi- 
circle ACB (b^^/b9) = 0. Therefore, in passing along 
any curve from a point h on the diameter to a point 3ST 

on the semicircle, the ratio { j :[^ \ passes through 

all possible values from 0 to <»; hence it follows that 

if soii^e value ft l/(2p+l) is chosen for then on 

each of the curves joining k and il a point will be 
found at which the expression 



becomes zero. The series of these points in the T^g 
region will be on a certain curve. The point' where the 
latter meets the curve T a= will be the point at 

which there holds the equation. 

VT.e/ VbV ST KhdJ 



and therefore ty the preceding formulas also the equation 
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Evidently from the manner in which one of these 
points is. ohtained it must he concluded that they form 
a certain dense curve. Thus '^(x^y) and 9(x,y) will 
not l^e single— ^''alued as is required in a real motion of 
the ;'^as, ThuSj_ i.f from the condJ.tion_s of the £rohlem_j^ 
it j.^ ££^^^"^1^ ^£ ££11^-^^*^^ IhSil, i'^£ pressure _of jth^ 
£i:-2i-^£2: Z£i:££i:^Z j P^Z'^^^i^l e>>:c e ed the limit 

^££l£;£^ jy^t^he"^Tnequal i jy T < lJX2p 4- ly , jthen st eadj 
mot ion i^_s, at^ least , n o t^alway s pos sTble , 



The author returns to the solution of the problem 



which was especially thought of in 
The flow o.f a gas hounded "by plane 
se:onrates and continues to flew in 
pressure is." cons idered. The prolyl 
out of a very large vessel and the 
gas flo\7 at a plate will he studied 



setting up this analysis, 
v/alls at v/hi.'ch the a'as 
a region of const an t 
ms of tjie flow of a fT^as 
pressure of an infinite 
in greater detail. 



Oonsider a 
•the form 



particular solution of equation (12) of 



= Zj^s in( 2nG + a^^) ' 



(14) 



where 



func t ion, the ordinary differential equation 



is a function only of J, To- determine this 



d r , 



>B dzn\ l-(?p+l)T 



dT J 



- — ~ (l— Tj n"z = Q (IE) 
t(i-^t) n 



is used, or, explicitly 



^2(1^^) t[i+ (p^i)T]l/fl^n^ Ci-(2P4.1)t]2j,== 0 (16) 

dT dT / u 



Setting 



zn = T^y , if n > 0 (l7) 
n 



yields, for the determination of y^^, the equation 
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2 , 

(l^T) [2n+l+(g-2n-l)T] -^ii + pn(2n+l)yn= 0 (l8) 



This is a hyper georaetr ic equation. Its integrals 
are of the fore 

where iC(T) denotes the series 

If it is desired to have an expression for that 
does not become infinite at the critical point ^ 0, 
in equation ( 14) , it is necessary to take the integral 
of equation (l6), v/hich remains finite for - 0, It, 
therefore is assumed that, "by making use of the notation 
of Gauss 

=^ ^^^n* ^n^ ^n + 1, t) (19) 
where a^j^ and are determined from the equations 

a^ + bj^ = 2n-P, a^bj^ = - pn(2n+ I) 

The question to be decided is v/hich of the problems 
of the above— indicated type may be solved with the aid of 
a fimction expressed by the formula 

^1/ = A + Be + ^^n^^n (20) 

v/here A, B, Bry are certain constants and is deter- 

mined by formulas (14), (l7), and (19). 

First consider the boundary conditions of the 
problems. Since the gas mass is to be bounded by stream— 
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lines, then along the bounding contour of the T^g 
region ^^ must assume certain constant values. If the 
part of the contour under consideration corresponds to 
a plane v/all, the angle 9 formed "by the velocity d i— 
rection with the x— axis should maintain a constant value 
so that th.is part of the "boundary will \>e a section of 
a straight line passing through the pole T ^ o. If the 
surface of the Jet is considered, there is along the sur- 
face p « constant, and therefore, "by the Bernoulli theorem, 
the velocity should likewise have a constant value v^^ 
But v2/2a = so that t likewise has a constant value 

T^^ It is clear that the part of the Tjoundary of the 
0 region corresponding to the jet will "be formed of an 
arc of a circle the center of which serves as the pole. 

The problem proposed of the motion of a gas mass is 
now compared with the corresponding problem of the flow 
of an incompressible liquid for the same boundary condi-- 
tions (the same disposition of the boundary walls, vel- 
ocities at infinity, and velocity at the jet boundaries). 
The latter problem is solved with the aid of the well— 
knoi/n Joukowsky method. By the use of this method the 
relation between the complex variables Ig v^/v + iQ = 

y '^o/'^ w cpj^ + i^^ is found where cp^ 

and are the velocity potential and the stream func- 

tion corresponding to the problem* It is assumed that 



is obtained and that this function can be expanded in a 
series of the form 




(21) 



w 




\ 

ie ' + Sk, 



n 




Then 



- A + B9 + Z3. 



'n 



sin (2ne + ajj) 



(22) 
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11 ll £s_ser_ted t ha t th_e £or r e s£ondi.n^ ££0_blem in the 
£^1® 2l i^i^ 11^^ 11 solved t h_e f ormula 

= A + Be + 2Bn f--^ sin(2n9 + an) (23) 

VoJ 7no 



where is the hypergeome tr ic series defined. by formula 



(19), y its particular value with T substituted for 

' n, 0 

and \ a certain constant. 



The correctness of this stateinent may, in part, be 
proved immediately. Thus, it is readily seen that for 
T = the right— hand sides of formulas (22) and (23) 

agree; hence if for t = = constant, then like- 

wise \|/ - constant. If, further, for any value 9 = 9^ 
the function aefiiied by formula (22) does not depend on 
this is true only if the condition sin( 2n6Q + a^^) 0 

is Ga.tisfied for every .n under the summation sign; but 
then the right side of formula (23) for the same 0 also 
will have a constant value. Thus the boundary conditions 
imposed on the function >V are satisfied. 

It is now noted here the the series formally 

satisfies equation (12)^ since it is the sum of its par- 
tial integrals. If^ £ow J.Jb J._s ^h^wn Jbhajt £or an;;^ JLJSLJ!fi 
Ih^ (23). £^nver^e_s an^^ ££i JL^-,-T£ JE^^^^ 1^ 

l^I^l llHiH 3.31.1^3. LSSl, then the function expressed by 
it actually will be the required stream function. If, 
moreover, it is shown that this series converges absolutely 
and uuiformly together with the series obtained by its 
term— by-term differentiation with respect to t and 9 
it will be justifiable to consider the latter series as 
expressions for the partial derivatives of the initial 
series. Then for a given making use of equations 

(7'), (80 and (11), cp, X, y, . will be found.. As regards 
cp, it is found from formulas (ll) which lead to the re- 
lation 



2sin( 2ne+an)dez <j^t( 1-,t) ^ 



.2n( cos2n0 + a_)(iTz^ iniimlZ (i^t) ^ 

2t 



-e-i 
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whence "by use of equation (15) and the following ones 
equation (24) 



= C + 3(1-t) 



2 J T 



•0 



dT 



(l-.T)"^^2Bn /Ciy^^IlL^ A-,JI!i^^ cos(2n9-H.a^) (24) 



may be easily obtained. 



The functions 



T y» 
1 + a 

n y^ 



which, .in what follows. 



in 



will be denoted by x^, play a very important part 

this problem since through them are expressed the char- 
acteristic constants of the various problems^ Certain 
properties of these functions and the methods for their 
computation and likev/ise the essential properties 



functions 'Zj, and y^ 



of the 

of interest here will be described 
in the following section. Only by becoming acquainted with 
all these properties is the possibility obtained of demon-^ 
strating the correctness of these statments that. remain to 
be proved with regard to the fundamental series for 
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PAET II 

CSaiAIiT PEOPSRSIES OP TH3 FUNCTIONS z^, , AFD x^. 
PEOOP Oy CONVEHGSIICE OP THE SERIES POH \|/ AUD cp, 

is that integral of the eq.uation 
n^Cl (2P + 1)t](1 t)""^*"^ (25) 

which does not hGcor.o infinite for T = 0. This integral 
is of the form Zn = ^^Jn where n > 0 and yn is the 
hypcrgconetr ic function 

= P(an, 2n + 1 , t) 

the parameters a^ and h^ hoing defined ty the formulas-. 

a^^ + hj^ 2n - P, si^'b^ = -pn(2n + l) 

It will ho shown first of all that does, riot 

possess any reaj. JJ^ots ^*wo£n tjie values 0^ 2p^~i 

of tho 5ca^J^Jlle_._ T • Assume the contrary to he the case and 
let T = a he the least positive root of the function Zj^. 
Since hecomes zero also for t = 0, then, hetwoen 

the values 0 and a, a quantity h should exist which 

d Z n 

serves as the root of the equation = 0, Thus the 

dr 

function under the differentiation sign on the left— hand 
side of eq,uation (25) will have the roots T = 0, 7 = 1), 
and consequently its derivative should possess a root 
T = c where 0 < c < h < a. In view of the fact that on 
the right-hand side of this equation the coefficient of 

cannot hecome zero for T < — it must he assumed 

23+1 

that Zyx(c) = 0; and hence the function z^ must have a 
root T = c < a. By the same reasoning it is concluded 



The famction Zj. 

T~T (1 - t) — - = 
dT dJ 



S4 
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that in the range of values of the variatles from 0 to 
a an infinite numl^er of roots of z^^ must be included. 

But then Zj^ could not he expressed hy a power series. 

Iron the proposition just proved, it follows that 
is an increasing function. Thus, since z^^ has no roots 

less than ^ , z- must always vary in the same sense. 

2P + 1 " 

Since it is pos-itive for very small T, the same sign will 
he maintained for all values of the variahle within the 
ah ove-ment i oned range. Thus z^ increases for values of 

T near zero, and hence will continue to increase until 

changes sign. It is noted that 

vhcre it is concluded that, £or 0 < T < ~~ — - , t he .f unc- 

Sp + 1 

for v/ithin these limits > 0. 

Ttirn now to the function yj^ , the holomorphic inte- 
gral of the equation 

y\ t(1 t) + y'nC(2n + l)(l t) + 0t3 + 3n(2n l)yj^ =: 0 (26) 

or its ec[uivalent 

J:.T^^+' (1 t) y' + pn(2n 4-1)7 (x ^ t)"" y - 0 (26 0 
dT ^ n 

Prom the theorem -just proved^ it is concluded that 
y^^ does not have any roots hetvreen the values of T with- 
in the ranjje ^J^^ sj.de rjd^ The same . may he proved likewise 
with regard to'the successive derivatives of this function. 
First of all, from equation (2S 0 it is concluded that 
what has heen started is true t'/ith regard to the function 
Thus, if there existed a root of this function the 
derivative ' 



0 < T < — i — . Thus ;-^„.> 0, 7' < 0. y"n > °' 

3P + 1 ^ ^ 

y'" <0 and the quantities , y' , y'' numerically 

decrease* 

■ Consider now the function s^. = ^n^^ ^ ^^^^ 

differential equation v;hich this equation satisfies can 

"be easily derived froTn (36) and is of the form 

r 

(1 « t)-^-.t^^--^' (1 - T)^^-P s» + T^^ (i ^ T)^^~P jOn - m)(2n + l) 

dT ^ . L 

+ n(m « p 1) = 0 (27) 

1 - tJ 

37 setting IT; = Pn it is seen that for n > 1 + ^• 

P 

the quantity within the "brackets maintains the plus sign 
whatever the value of T*. If n is an integer this is 
true for all n > 2; for n = 1 it \rill have the minus 
sign. It is assumed that n > 1 + . 

3,y setting m = |j*3n for a suitable choice of [j, , 
there is in the hracliets a negative quantity for all values 
of T if the expression is negative or zero for the larg- 
est adiiiissitlo value of t; namely, — . ^or this it 

2p + 1 

is sufficient to choose IJ* so that it satisfies the equa- 
t ioia 

(1 - li)(2n + 1) + -li- (i^pn - p - l) = 0 

2P 



or, after reduction, 

p,^n - 2\x (^2n + I + J--^ + 2(2n + l) = 0 (28) 



*lfithin the liiriits, of cout^B^, of 0 and ^ 

2p + 1 

this- uust "be kept in mind t hr-^5\:^hout the f ollov;ing discus- 
sion. 
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dT ^ 



would likev/ise have a root. But this is impossiljle "because 
the second term on the left-hand side of" equation (26 0 
cannot hecome zero within the range of variation of T. 
By differentiating m times equation (26), there is oh- 
tained 

y^^(m+2)^^ -7)4. 7./^^'^^^[(2n 4- m + ' 1 ) (l - T ) - (m ^ $)t] 
+ [0n(2n + 1) - n(2n - 3) - m^]^^^^"'^ = 0 

or 

-fi^T (1 - t) 5- + [Pn(2A + 1} m(2n - p) 

dT n ^ 

- ^s^^n+m ^ ^ ^)m-^P-i ^^(m) ^ q 

whence r eas oning as "before, it is concluded that yn^^"*""^^ 
cannot have roots within the range of variation of T if 

y (nl) does not have roots within this range. By setting 
n 

m = 1, 2, 3 the correctness of this statement is 

pr oved. 

Setting T = 0 , gives 

(0) = 1, 7\{0) = ^pn, y" (0) = rpn(2n + l) 2n + p ^2 

^ ^ 2n + 2 



(0) . - tl^JtJljzJ^.l±^J:i^r.i2n 1) ^ 4n ^ 2P - 4] 

^ 2n +'2 2n + 3 



Since P is approximately 2.5 and n is a positive 
number, the signs of the .above quantities alternate. The 
same will he true for any t satisfying the inequetlit ies 
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The maximnni value of |x is olitained for very large 
n' (n = 00 ); in this case |x = 2 . In general |a» < 2 ; 
thus, for n = 2, P = 2.5, |x = 1.083; for n = 3, 
Ui = 1.181, 

With the a.lDOVG choice of the quantity m the func- 
tion cannot have a root different from zero v/ithin 
the ran/^e of variation of t. This can "be proved "by the 
same reasoning already more than once applied. 

It is noted that ^^^(0) - 1 and from equation^ (27) 
s*j^(0) = m - pn is ohtained. For m s= j^-pn, s'j^(O) > 0 and 
therefore the function 

y^d - t)-^^*^ (29) 

increases with the variable and v;ill exceed unity. If 

m = however, s» (0) = 0, TDut s" (0) = -3n 

^ 2n + 2 

and this magnitude for n > 1 + — is a negative quantity. 

P 

IT or this reason ^'^^'^^ likewise, as a decreasing function, 

should for t > 0 be less than zero; hence it is concluded 
that 

yji - T)-P^ (30) 

for n > 1 + i is a decreasing function and represents a 
P 

proper fraction. 

Por n < 1 + - the quantity within the "brackets in 

equation (27) will "be less than 0 for m ^ pur s*ji(0) 0, 
s"n(°) > 0» stnd therefore ' s ' ( t) > 0 and ' 

y^(i - T)"-^^ (29 0 

will he an increasing function. 

The smaller root, denoted simply by |J» , of eo^uation 
(28) yill be less than 1, It can' be readily shown that if 

m = [iPm and 0 < T < — ^ — the coefficient of s^ 

2p + 1 ^ . 
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in eq^uation (27) wil3. again "be of con5ta.nt sign, the 
latter being positive. Therefore s'j^Ct) > 0 is obtained 
since s^i(O) > ^5 hence it is concluded that the 

f unct ion 

yjl _ T)-*^^^ • (30.) 

decreases with increasing T. Thus, for example, for 
n =^ 1 , M» =-0.93 (P is taken equal to 2,5) , and there- 
fore with increasing yi(l t)"^^ , yi(l - t) will 
"be a decreasing function. 

The above-mentioned properties of y^^ give limiting 
functions within which y^i is inoltided; namely, 

for n > 1 + (1 
P 

for n < 1 + (1 



,Pn . .\^&n 

t) > > (1 - -r) 



^) < yji < (1 - t) 



where |Ji» is determined by equation (38) and is equal to 
the smaller of its roots. 

It may be noted further that the function 

ynd - T) 



increases a fortiori, For n < 1 + ^ this is evident; 

for n > 1 + -i, on the basis of what has been said above, 

it may be considered as the product of two increasing fxxnc- 
tions y^il t)--^^^ and. (l - t)-^^-^^^^^. Hence if the 
greatest value v/hich T assumes in the given problem is 
denoted by and the corresponding value of yo by 

jr^ ^, . this will give the inequality 

' ^sPn 
T::^Pn < 

or ' 
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Ji_ < ^-=-11'^" (31) 

T y ' -n 

Last to "be considered is the function x^. = 1 + 

. • . ^ yn 

on which , as has iDeon shown in the for e going , depends the 
computation of the very important constants of interest in 
the var.ious problems. The differential equation which the 
function satisfies is first set up. It is obtained 

.froLi. the hypergeometr ic equation (26) by setting 

T 

xn- 



= e 



Thus, it is found that 
x»^t(1 - t) + lix^^d - "T") + x^PT" - ^-[1 - (23 + 1)t] = 0 . (33) 



This equation , together with the condition 2:^^(0) s I, 

fully determines the function x^^. It has "been shown 
already that the function Xn for a change in the variable 
within the limits under consideration remains always 
greatqr than zero. It v/ill be shown tha t it de creases with 
increase in T. P or this purpose equation (32T was differ- 
entiated; there was obtained 

x"j^t(1 - t) + x»n[2^(l - '^)^n + (p - 1)t f 1 - t] 
= nx:i^ - ^ n(2P + l) 



Substituting in the brackets for x- its value 1 + -1 

a a ^'^^ 

and multiplying the equation by T^^(l - ) ^n » reduces 
it to the form . 
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T y * 

But is less than 1, since x„ - I = - — - is a neg- 

^ ■ n 

ativG niagnitude on the "basis of vhat has "been said with 
regarc^. to the signs of the function and its deriva- 

tives* Hence, the right— hand side of the obtained equa-^ 
tioa is a n-egative quantity of constant sign. If x'j^i 

equal to — ^ for T = 0, had a root within the range of 

, 2n+i , -3+1 3 

variation of t, then -f^j (1 -t) y- x' would 

also tecome zero for a value of T less than this root, 
a result v;hich is impossible. But x^n» everywhere 
finite, as can readily "be shown, cannot change sign except 
"by passing through a root. Thus x'j^ remains less than 
zero and therefore decreases. 

The next step is to seek to olDtain functions that 
limit the value of x^. For this purpose the following 
theorem T^e proved: If^ on substituting in the e g^ua - 

t i on determ ining Xj^ , a holomorphic f unct i on , ther e 
is obtained on' th_e_ lef £ ^s^ide^ .£0.siti_v_e Y-?il^-?...?X 
sign, then kj^ > x^i; the inequality sign will be reversed 
if the result of the s ubstitution is less than z ero , P or 
T - 0 , may be equal to 1, From the assumed inequality 

jjT(1 T ) + nkj./(l - t) + k^pT ^ n[l - (2p + 1)t] ^ 0 (33) 

subtract equat ion' (32) , which leaves 
J,*- x\)r (1 -:t) + (k^ ^ Xn)[3T + n(l - r)ik^ + x^)]^ 0 (S3 0 



By setting 



T 

k.-i 

ndT 



= 1 + l i^S; ^ > 0 

n tn 



then, on substituting in the brackets for k^^ and x^^^ 

their values in terms of In and y^ ^nd multiplying by 

—3 

the positive factor yn^n""^ (l - t) , there is oTd- 

tained 
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Integrating this 
yields 

(kn 

whence the rectuired inequality is oTjtained 

^ ^a 

It \7ill now "be shown how the function niay iDe com- 

puted to any degree of accuracy hy transforming it into a 
continued fraction. In differential equation (32) the new 

T 

independent varial^le s defined loy s ^ T~""f" substi- 
tuted; when T varies from 0 to 2p-^+~» s varies from 
0 to ~. The differential equation for ' will bocome 

x'^ s(l H- s) + Xj^Ps + Xix^^ - n(l - 2Ps) 0 (S4) 

From equation (34) is found x^CO) = 1, x'n(O) - "^P • 
any function satisfies inequality (33), then on sub- 

stituting in (34) there will be obtained 

k^j^ s(l + s) + kj^Ps + nkj^^ - n(l - 2Ps) < 0 (35) 

whence follows as before the relation k^ > Xj^ , for 
kn(0) ^1. 

Equation (34) together with inequality (35) will be 
writtenas. 

x^^ s(l + s) + Xj^N + nx^^ - n(l - 2ps) > 0, (36) 

which is to be under st ood as follo\\rs: If, on substituting 
any finite function within the range of variation of s 
and equal to 1 for s = 0, the result is zero on the left- 
hand side of relation (36), then this function is the exa'ct 
expression for x^; if , as a result of the substitution- 
a positive quantity is obtained, the substituted function 



inequality within the limits 0 .to T 

^ sn . , > ^ 

- xn)T (1 - t) yn^n < 0 
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is always greater than x^; in the contrary case the sign 
of the ineq.Tiality is reversed. 



(37) 



where c^^^ is a new function to be determined. On reduc- 
ing and changing sign, there, results from (36) 



2 



3(1 + s) + c^""^ (2n + 2 ^ Ps) - (2n + l)c^^^ s ^ nP .+ 0 + 1^6 (37O 

where it v/as necessary to reverse the inequality sign. 
The meaning of the relation is as follows: If, after siib- 
stituting in the left-hand side any function in place of 
c^^^^ the result is a negative (quantity; then replacing 
c*^^^ in formula (37) hy this value there is ohtainedji^n 
upper .liniit of t_he 5.1.^5^51 . .^ji - that is , a function 
greater than 

Further is set > v 

(n) 



in) ^ ^ 

1 &'"'s 



c =s -5— -y— (38) 



where c^^""^ = c^^'Ho) = ^ - ^^-i-i , and the function h 

° 2 2n + 2 

satisfies the relation 

h^""^' s(l + s) + a^^'^CSn + 3 + 0 + Ds] - h^""^^ s(2n + 2) 

- (2n + Dco^^^) - ? 5 0 (38') 

whence is obtained 

^0 - ^ 2 2a + 2 2n + 3 



Next, setting successively 
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yields, for the determinat ion of Cj^ , , Cg, dg,.., 
the relations: 

c/'')'s(l-|-.s)4-r/"»(2w4-4— ^s)— (2w-|-3)c/">-s- 

— ^„"''(->«-f2)-f(i-hi§(). 
(), '"''.sfi +s)-f (y,("»| iH-f-a+dS-fijsj— (2m-|-4)(>, 
•_c,.;'''(2«-h:))-M0, 
1 ^s) 4-f ^(")(o„ ^6— ;3s) -(2m -f 5)c;"''-.<»- 
-,V.r(2«+4)+fl4-l5 0, 
,,/'')'.s(l-j-s)-f,//')[2y^-|-74-(|i-f-l)s]_(2n-l-6)'V'''-s- 

-<',/''(2«4-'»)-3§o. 

HoreoYer 

e(").^^_L02i+i 42^-1 ,42^+1 

^^ - 2 + - -t.j,,4_4 ' ',,0-0 + ^^2^ - *'2ii4:5 ' 

With the aid of these formulas it is not difficult 

to set up equations for the determination of 0^"^)^, 

and formulas for c^*) . d^^^ The 
■» m, 0 • m, o 

latter are of the form 

op^ after reduction: 

"'•" 2 l-HT^ '(2» + 2/H-f l)(2H-f2m-|-2) ' 

(40) 

^c; _ ^ ...g, n 0>H-l)(2n+m+l) 

2 ^"'^ ''f2w-f 2/H-f 2)(2>/-f-2«/+.l) 
The equations for determining Cj, and are 
the following: 

c<"^,s(l+s)+c(''',„(2n+2»«+2-p«)-(2«+2w-|-i;c<''>\,s-- 

^"V(l+s)+d('^,,[2n+2;«+:H(|3+l)>«']-(2H-2/«+2)t^<''^,'- 
-c''Vo(2«4-2m-f 1)-13§0. 
All these relations can oe readily verified by the 
method of passing from m to m + 1. 
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Ihe sign < in the last of the foregoing relations 
holds until the index m exceeds a certain limiting 
value: namely, while d^l^^is positive; if c^n;^^^<o 

then in the relation for S^^^j^ the sign < must re- 
placed "by >. This is hecause, among the simplifications 
which were made in transforming the ahoVe relatiQft, there 
occurred division by ^^^^u,o' 

By collecting the results, x^^ is finally expressed 
"by the formula: 

, Ps 
= 1 - 



1 - c 0^ 



1 - ^^""^s, (41)' 

""~Tn7"~ s 
1 . 0 

— --(ny— - 

1 - ^ a.o^..-- 



where c^^^ and ^^^^ „ „ are expressed hy formulas 

(40). 

(n ) 

How consider the magnitude of the ctuantities c m , 0 ' 
^(11)^ 0* difficult to see that thoy are always 

contained hetween P/2 and -X. The first of these is 
ohtainod for n very large (n = ») ; the second, in general, 
differs little from -l/4 and is ohtaincd from the minimum 
of the expression 

(n) ^ I M. l)k(2n .+_kl__ 

° ^-^-^ 2 2(2n + 2k - 1) (n + k) 

as a function of k. This minimum occurs cither for 

k = Eny4n'* - 1 + n(2n - l) , or k = ^/ An^ - 1 ' - 

+ n(2n - 1) + 1; these values for integral n are equal 

to 

4n^ _ n - 1 and 4n'^ - n 

2 B 4* 1 ( n } 

The coefficients of —-g in the formula for c k-i,o 

correspondingly receive the values 
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and 



2 

2x1-1 

2 2 

(4n^- - n - l)(4n - 1 + n) ^ ISn" - 1' 



(^.-:rK) 



the second of which is larger than the first, Thus the 
mininuni value of c^'^'j^.^j o 

3 ^ £|_t_i 

... ^ 



4 1 



^ 1.6n^ ^.1 ^ 



It can further easily te shown that the continued 
fraction (41) is always convorgent. The contrary could 
"be the case only if the expression 

c(-)oS 



1 - a'^^u 



1 - c"^^^ 



0' 

TnT 



approached unity. But even in the leas t ■ f a v orah le case, 
for n = 00, this quantity T^ecomes 

!!!__ 

1 - i$s 
2 



and its raaximum value, obtained for the maximum value s 
is 
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As regards tho signs of c and h (for simplicity 
in writing, tho indices are omitted) those for 
which m is equal t o 0 , 1 , 2 • , ^ up t o a certain limiting 
value will he positive, all the remaining ones negative, 
provided that n has a finite value. The limiting value 
of 1 + m is olDtained from the inequality 

i. (2P + 1) ^l^^l < 0 

2 (2j^ + 2k - 1) (2n + 3k) 



or , on reducing , 

k^ + (2n + $)k - Pn(2n - l) > 0 

This inequality is satisfied as soon as m + 1 = k 
exceeds the larger of the roots of the equation 

+ (2n + p)a - 3n(2n - l) = 0 

The limiting value of n will therefore "be expressed hy 
the f.ormula 

jn « B (^n - £ + yn^ (2p + l) + £f 1 
I 2 4 j 

If only integral values of n are considered, the 
following limiting m and ^^^^m^o "b© o*btained: 

n = 1, m = 0, c^'\ « 

n = 2, n = 1, c^^^ ^ =: Ji- 

I'O 28 



^3 ) • • 17 P assumed equal to 2,5 
n = 3,. m=3, C'a.o" 

n = 4, m = 4, c^^\ = ^ — (42) 

3X4X7 

m inci-easing with n as req.uired. The quantity ^^^^m o 
likev;ise becones negative, "but only for large values of 
m, Por ^^^^ni.o *° ^® less than zero, the inequality 
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£ _ (2P + 1) ^I2n + kj < 0 

•2 (2n + 2k + l)(2n + 2k) 

must "he satisfied, or , 

+ (2n - P)k - Pn(2n + 1) > 0 

The limiting value of m denoted "by nii is expressed 
"by the formula 



= bI _ n + £ + yn2(2P + 1) + 

I 2 , 4 j 



But if p > q. and p = S(p) + e, q = E(q) + ^ 
where 9 and ■b are proper fractions, then p - q =.E(p) )-.3!(q^) 
+ 9 - •6 ; thus 

Il(p - a) = E(p) - , or S(p) - - 1 

Therefore, comparing the obtained values of and 
the limiting m, it will "be found that m, will te ©q.ual 
to the limiting m for the coefficient cvn)^^^^^ plus 

B(3), or plu.s EO) + 1. Thus 



n 


1. 


mi 


= 2. 




0 


n = 


2, 


mj 


= 4, 




_ -L. 
"* 28 








5. ■ 


0 5.0 - 


_ -L 


n ss 


3, 


mj 












76 


a - ■ 


4, 


.111 


= 7, 


0 7.0 - 


3_ 

100 



P = 2.5 (43) 



All c^^^jji 0 starting from that which corresponds to 
the limiting value of m like all ^i^^^m.o m ^ mi 

1 — 0 ffii, 0^ 

are negative (quantities and the iiemai^nde.r 7— y 

1-c ^^^nii+i .oS 

of the continued fraction is expressed in the usual form 
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1 "f a 

1 + Tb its numerical value "being contained 

1 + c * . • 

tetvieen •! + a and 1 +~a * 

ITT -b 

Of the functions Xn the one that is particularly 
simple is x., ♦which for the assumed value of P is ex- 
pressed as a'fraction of two polynomials of -the third 
degree • Por n = 1 

e(') - ^ S^^^, - 0 (by (43)) 

2,0 14 s , 0 

X, 1 ^ ._a^/2_ 



1+6/4 

1-7 s/lO 

1 + 7 s/a 0 ^ 
1 - 3 s/a 8 



1 + 5S/l4 

or, after reducing, 

X ' - 32 64s - 14s^ - 8s^ ^^^j 
V (4 .+ s) (s^ + 2s + 8) 

With this formula y^ is readily found* Tor this 
pufposQ the previous variable t a • — is substituted 

1+8 

which gives 

y'l 32 - 160t + SIOt^ - 84T^ . . 

Xi = 1 +t1~A = _-- . (44 «) 

yi 32 - SOT + 70t2 _ 2it3 

vhence, s1i,pce yi(0) = 1, there is obtained 

32yi = (4 - St) (8 - 14T + 7T^ ) 
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The preceding simplification results only from the rounded 
values: 3 = 2.5, Y = 1*40. The more accurate value of 
Y for air is 1.4025 and 3 is somewhat loss than 2^5. 
In this case all Xn Q^re expressed Tdjt infinite continued 
fra.ctions, 

^ To employ in the applications the exact formulas for 
Xj^ appears impos s it) le , since this would offer very great 

difficulties which have not heen overcome* However, "by 
usiac even the simplest proper fractions, Xn is obtained 
with sufficient accuracy. Entirely satisfactory results 
are ohtained even in the case limited to the' third proper 
fraction and x^^ is expressed "by 

= 1 - Ps 

1 - ^Jj'lol (45.') 

or, after reduction and suId s t itut i on of the values of the 
coefficients 

^ Ps^(2n + Ll«il (45) 

^ (2n + 2)[2n 3 - (Z^n - 3 - -2)3] 

The error for such computation of Xn is greater the 
greater the value of s. The magnitude of. this error now. 
is estimated, cons ider ing only the_ integral va lu es of n. 

With the exact value of Xj^ for n = 1, a direct 

comparison of the results of the computation of this func- 
tion may he carried out. hy formulas (44) and (45). The 

computation will "be made for maximum s = or s =5 0^2, 

2p 

assuming as "before 3 = 2,5, The exact value of xi will 
■be 

7 7 S 

= 0.52538S3 

1477 



Prom formula (45) Is obtained 
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Xi ^ 0.525510 

The difference is approxiiaately 0.00012'. 

Tor other values of n in estiiiiating the errvr , it 
is necessary to proceed otherwise. It is noted^fltst of 
all,. that for all n > 2 formula (45*) gives a function 
greater than x^i the contrary is true for n = 2, In 
order to show this, turn to equation (33'), deter laining 
Shis equation may be written as follows: 

h^""^ 's(X + s) + 6^''^[2n + 3 + (P, + 1) s] - a^''^^s(2Ti + 2).. 

- (2n + 3)b^^\ 0 

• Substituting in the left-hand side for give 

b^^^o^Ce + 1 - (2n + 2)^^^^ J 

or, after substituting in the brackets for S^^^o -its 
value by formulas (40) 

)^sri ^ Pn ^ l2P^^l)l2nJ:^l) 1 ^ ^(n) ^ ^ p^^^^j ^.^MjL 
L 2n+3j °L .2n+; 

Por n > 2 this value will be less than 0. But from this, 
as has been said, it must be concluded that on replacing 
h^^^ in the formula for by a trial value, a function 

greater is obtained. On the contrary, for n = 2 in 

the brackets, the. quantity 

2 - B^^^t^^ > 0 
7 

which shows the correctness of the reversed inequality 
(xg is greater than the value that would be obtained on 
substituting S^^^q in the formula for x^) for 

the value n = 2 . 

As regards . the lower limit of the functron xirill 
be the following proper fraction. Thus 
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3 (S) 

Ps Ei_£ 0 < xs < 1 - Pb 



,2 



Ps^c 

1 - c(^), 



1 *o 



By computing the values' of the limiting functions for the 
case of the greatest difference in their values, v/hen 
1 

s = -— a 0,2, there is obtained 
2p 

0.47C34 < XsCO.S) < 0,47037 

(n) 

In the case n > 2, assume in the equation for 
/_N ;i(n) 

s ~-- — .a, where k is a constant to "he determined, 
1 - ks 

and k is chosen so that the result of the sulDst itut ion 
is 2^<5atGr than zero for 0 < s < g^. This requirement 

leads to the inequality 

k(2n + 4 - 3s) - k®s(2n + 3) - • (3n + 4)c^^\ „ = 0 



The smaller root on the left— hand side of this xn- 
eq.uality is expressed hy • 



^ + 2 " - VCn + 2) a - 2ps) + 4s (n + l) (2p + l) + 

2 4 



(2n + 3)s 



Its maxiimm value corresponds to s = -1- = 0,2 and is 

2^ 

equal to' 



n + 1.75 - y2ilB_±al + 
2n.+ 3 
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If k is equal to this value., the ahove inequality will 
"be sati-siied. 

It is nov/ possible to indicate the limits within 
which is included. The upper limit is expressed 

fornula (45'); the lower is ohtained oj substituting in 

this fornula — .-Z-O ^ h q. After reducing, finally 
. 1 - ks 

Ps^On - p - 1) 

1 - Ps > Xji >' 1 - Ps 

^ 2pn - p _ 2 V ' 
(Sn + 2) II s ) 

2n + 3 ^ • . . 

3s^(3n - p - 1) ^ ^ 
_ (47) 

, , 2Pn - p - 2 ,\ 

(2n + 2) ( 1 s - ks^) 

^ 2n + 3 / ' 

where \ = 5 = i - i££_i_iL(£^_l_il and k 

1 - ks ' ° 2 (2n + 2)(2n +3)' 

is determined hy (4o). 

She iiumerical values are given for s = 0.2 of the 
limiting ff.nctions f or x^ for n = 3, 4, 5, and 6: 

0.4348 > x^(0.2) > 0.4343 

0.4095 > X4(0.2) > 0,4073 

0.3905 > X5(0.2) > 0.3872 

0.3755 > Xg(0,2) > 0.5704 

It is thus seen that the error increases, or more ac-curately 

may increase with- n hut nevertheless is very small for 

small vai^iQs of the latter. For somewhat large values cf 

n the li":its of error v/iden. Thus • 

n = 12, 0.326 > x^^(0.2) > 0.308 
n = 24, 0.293 > X24(0.2) > 0.251 
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This unfavora'ble circumstance is offset, however, to some 
extent hy the fact that the functions x^^ with largo n 
enter the inore removed terms of the series and the coeffi- 
cients of these terms are relatively small. 

The limiting functions for x^ also will he given 
with large n. These functions will he useful in computing 
the limits within which the remainder tern of the series 
for the gas Jet prohlem is included. Again, in the differ- 
ential equation defining xn 

x\q(1^ s) + x^Ps + nx^^ n(l - 2Ps) ^ 0 



Substitute on the left— hand side the expression 

y W"3|3 G + 3 us^ , and choose the function u so that the 
result of the suhstitution is greater than zero. Then, 
hy the theorem proved ahove, 

- / ■ - 

V 1 - 2ps + 2us > Xj^ 

This, suhstitution gives on the left-hand side of the equa- 
tion for x^ . the expression 



u^sd + s) + 2uLl + (l + P)s] - 3(1 + 23) + 2nuyi - 2ps + 2us^ 
k 5= • ■ — s-" (Uy) 

/l.~Tgs + 2us^ 
which, as can'^Be seen, Mill he greater than siero if 



2n^ 2" n 



Thus, for this value of u, ^ 1 - 23s + 2us^ is a decreas- 
ing function of s; the product, however, of this root hy 

2nu for the maximum value of the variable s = — is 

23 

equal to (l + 23)3, and ther e'f or e " t he quantity k remains 
positive. On the other hand, it becomes negative, whatever 
the value of s , if 

2n + 2 
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since the numerator in the expression for k in this case 
is equal to zero for s - 0 and, as a decreasing function, 
will be less than zero for s > 0. Thus 



2n® ^ , n H- 1 

By raising somev;hat the upper limit of the function 
Xn, the first part of the douTDle inectuality also can be 
transformed into 

n * 

and therefore can be expressed by the f oriaula 

X, - vT^B * x.y5!iiZil (49) 

n 

where is a proper fraction. 

It can be easily shown that x^ for the same value 
of the variable decreases with increasing n» This is 
clear from the eq.ua t ion for Xn+j^: 

'n+ins(l + s) + x^^^&s + (n + rsi)x^^^j^ ^ (n + m) (1 - 2ps) := 0 

.Substituting on the left Xj^ for Xj^+jq , there is obtained 
on the basis of the eq.uation for x^ 

EiCx^ii- (1 23s)] 

a magnitude greater than zero due to inequality (48) , and 
therefore it is concluded that 

> Xj^+jj^ (50) 

whatever the positive number m. 

iTow v/ith the properties of the function Xn that 

are of importance for what follows, two inequalities 
which the functions y^^ must satisfy will be noted 



X 
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further. She first of these will "be derived in the follo'.r- 
ing nr^naer : Set up the differential equation determining 



y'n 

^n = 

^ n 



(r\\ + r\^^)j(l - t) + [(2n + l) (1 - t) + PTjTij^ + 3n(2n + l) = 0 



2Pn 

Also set ^ 1^ ~" ^ » which gives 

1 T 



n u * 



« Pn(2n + 1) + -^^^I(23n p ^ l) = 0 (51) 

1-T 

To this equation, as 'can Tdg easily seen, the theorem 
proved for equation (36) for the function xn is appli- 
catle. If on suh s t itut ing for any holomorphic 

function there is ohtained on the left-hand side an 
expression greater than ^ero, the substituted expression 
will he greater than Ij^. If 2pn is substituted this result 
in fact is ohtainod, and therefore ln< 23n. 

Substitute, further, in the equation defining I n+m 
in place of In+m function |n* Il'^e result of the 

substitution, on the basis of equation (51), v;hich is 
satisfied by. roducos to 

1 - T 1 - T 



This expression is negative for any T, since it gives a 
result less than zero on substituting for 1^ greater 

magnitude 2Pn.. Hence tn^ tn^mt 



35i 1 - T y^+m 1 - T 



Integrating this inequality from T to Tq, and passing 
from logarithms to numhers, gives 



46 



KACA Ik Wo. 106? 



where o» "'^5^ assumed notation, denotes yn^'^o^* 
Thus the fuLction ^^n^ ( ™) docraase6.,xit&-J:S" 

y^^^V^X - Tq/ 

The second inequality which it was proposed to derive 
follows from relation (50). Trom the latter is -oT^tained 

y'n ^ y'n+m 
yn yn+m 

Integrating within the limits s and Sq, gives the 
Eesult. on passing from logarithms to numlDers: 

J^n^ < .Ill±5^ (53) 
yn,o yn+m,o 

yn 

Thus, the ratio — increases with increasing 

yn , 0 

It is necessary to proceed to the pr oof_o£,.the .con- 
V£r£;enca of the series giving the solution of the gas jet 
prohiems. In explaining the general method of solution 
of this type of proljlem (see pt , I), the following for- 
mulas for expressing the stream function and the velocity 
potential were arrived at; 

\l/ = B9 + y (~ ^ Jl&^ sin (2nG + a^^) 
^ ^ ^ ^ yn,o 

(1 - T)"^ y (-^-.7"' -^-Xr, COS (2ne + a^) 



-I) 



(52) 
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These are formulas (2S) and (24) of part I. The indices 
n eateriiit^' them increase as the terms of an arithmetic 
pr o;.;r e s s i on ♦ It v;ill "bo shown the preceding series are 
alDsoliitelv and uniformly convergent for any T < Tq if 
this is true of the series 

= B6 + ^Bn (^^- V sin (2n0 + a^) 
9i = B ^ ly - (^^ y cos (2ne + an) 

expressing the stream function and velocity potential for 
the corresponding problem in the ca,se of incompressible 
licLuids. The series and ^>l will evidently be 

lim Bn, 

absolutely convergent if < 1, where n and 

n = cc Bn. 

are tvro successive valtics of n. It will be assumed that 
this cor.dition is satisfied. On the basis of relation. (31) 
it can then be stated tha.t the terms of the series \i/ are 
correspondingly less than the terms of the series 



y [BnJ ( 



fjjLZjjjiy (54) 



which is an £;^^^l}itQly convergent series for T < Tq if 
Tf. < " ; for in that case 

4.-^(1 _ T)^^ = [1 - (2p + 1)t](1 - T)^^^' > 0 
dT 

' t(i t)^^ 

and therefore ^ ~ — - a proper fraction. 



The remainder term of the series \}/ 

-^n = y Bn Q-Y --a- sin (2ne + a^) 

n=n 
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is numerically less than the term R'j^ of the series (54). 

ns2 n " 

tut R ^jnt approaches zei'o with increasing n no matter 

in what manner T is less than • Prom this it is 
concluded that the series \|/ is iin if or mljr onvs.r £g.Gi • 

Since the series entering the function cp differs 
from the one Just considered hy having cosines instead of 
sines and the successive terins multiplied hy a series of 
deer eo^G in^i positive q^uantities, the theorems just proved 
likev/ise hold for the series cp. Turthermore, it can iDe 
easily seen that the same properties are possessed hy the 
derivatives of the functions (p and with respect to 

6 and therefore also their derivatives with respect to 
T, since the latter are connected v^ith the former hy the 
linear relations (ll) of part I. A consequence of these 
theorens is the continuity of the functions cp and \|/ 
and their derivatives within the rangfe of the variables 
T, 6 under consideration* (See vol I, p. SIO of refer- 
ence 7«) 

It will be shown, finally, that ^s . , ...T ^ „appr oaches t^h e 
ii^iiill£..v§ii2ie... .T Qj.,..iiie_^ .^aPEr oach limits 

,9lLS., -.2£ .^b ta ine d on . s ub s t t u t i ng 

"^P • this purpose, consider the sums a 'auid 

CT' "'of p terms after the nth in the series cp and \|/» 
Let ill, ng n be successive values of n; then denot- 

ing the ££aj£iiioil^ 



JL(^LZ™^>) y^ / 1 - T X -8pn 

To\l ^"^T^y ' yn,o^ 1 - ^0/ 



correspondingly by i and , these sums can be trans- 
formed into 





^nt^^n sin (2n9 + an) 
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The expressions 

cos (2ne + a^^) , |^ sin (2n6 + a^) 

will now "be denoted "by '^x^ and -u'ji. 

The series Sun and Eii'n uniforinly converge for 
any I less than unity, since they coincide with the 
series ente'rir.g the functions cp^^ and It is assusned 

that they converge also for | = 1; then from a knov/n 
theorem in .analysis their value for | = 1 is the limit 
which they approach as t approaches 1 (and hence 

But. in this case n can he given an increasingly large 
value so that thq sums ■ ' * 

^ni» ^ni^+ ^ns* ^n^ + ^ng ^ns » . . . ^n^ + u^^ + . . . + 

are included hetv/een any values e and 6 as small as is 
desired, xirhether the quantity | is less than or equal to 
unity. And, since on the "basis of the properties investi- 
gated in this section of the functions £^rid y-^, the 

quantities x^^ entering the expressions a and .a^ 

are greater than zero and decrease with increasing n, 
then loy the theorem of Ahel, cr is included "between the 
limits 

- (1 " t) Srij^^ and - (l t) cr\^^ 

Por the same reasons cr ' is. included "between other 
arbitrarily small numerical limits, and the propositioia is 
thus pr oved. 

As a result of all the properties which have been 
demonstrated of the series cp and \|; the conclusion is 
arrived at, which was the ohject of the investigations: 
namely, that the ..formulas, ^h^a ined....ar e_an_aciu 
i^G. S^^- ^ISM pr oTb lems , 
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PART III 

THE FLOW OF A GAS FROM AN INFIlItTELY WIDE VESSEL 



The method descrihed p-ill be applied to the prohlem 
of the fl077 of a gas from an infinite vessel ^ith plane 
Trails, the simplest case being considered — that is, 
where one wall is a continuation of the other. 

Consider an incompressible liquid flooring out of 
such a vessel (fig. 4); AB and A*B' are the traces of 

the walls of the vessel; OX is the trace of its plane of- 
symmetry; BCC*B' is the escaping jet» If the quantity 
flowing out per second is denoted by Q, the velocity 
potential and the stream function, respectively, by cp^ 
and considering \i/ = 0 on OX, then in the region 

of flow cp varies from *- oo to '+•00 and \l/ from 

0 0 
^ ji to + The complex variable 

3 3 

w =5 9i + i\|/i 

will be connected with another complex variable u (ref- 
erence 2) through the relation 

^ = S ~T 

The region of variation of w then corresponds to the 
upper half plane of the region of u. In addition, the 
logarithm of the ratio of velocities Vq/v at the jet 
surface and at the point of the fluid considered will be 
denoted by ^, and the angle of the velocity direction 
with the X axis by 0. The problem is then solved if 



d + i9 = i arcsin - 

u 

Q 

Tor, on the boundary ABC >|/ = - ^, c;* varies from -oo 
to +00 ancl "u passes through, the negative part of the 

real axis from Oto~cr; e = -,<»>t)>Oif 

2 
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0 > n > -c, ^ - 0, V = Vq, ~ > Q > 0 for -c>ii> 



On the ■boundary A'B'C* ^ ~ g' ^ varies within the 

same limits; e=--, oo>i>>0 for 0<u<c; t>=0, 

v = VQ,-^<e<0 if c<u<co. The point u = 0 

thus corresponds to the infinitely distant point of the 
vessel r/here ^ = 00 and the velocity "becomes zero; u = 00 
gives the part of the jet at infinity. Finally, for 
9-0, u is purely imaginary and ^\} - 0 and the center 
line of flo^j coinciding with the X axis is obtained, 

From the preceding formulas is found 

V7 = — — Ig ( 1 sin 

TT \ i 

or, if + iS is denoted hy cr, 

T7 = ~ Ig 1 sin - = Ig = (cr ^ Ig 2) 

TT ^ in 2 rr 

- a Ig (1 ^ 0-^^)* (55) 
TT 

It is noted that cp^ = 0 . at the points where the Jet 
separat es from the walls; cr in these cases has the 

value +i|. 

By expanding the logarithm in formula (55) in a 
•series, there is obtained: 

' n i) 



17 Zi 1 = ig 2 - ( ^ + ie.) + ) S (cos 2n9 - i sin 2ne) 

n = x 

whence 



*In the case of a vessel with the walls meeting at 

an angle it is necessary to replace in this formula a 

by ^ and in the succeeding relations by — , S by 

9 ^ ^ 

^. The angle between the walls in this case will be equal 

q 

1 0 qTT. 
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„ v|/i = -9 - \ 5 sin 2ne 

Sti"bst itut ing in this forinula' the variable T deter- 

mined by the equation — = = e results in the 

^0 V 0 

required expression: 



2n9 

0 



Since this series is aTDsoXutely convergent therefore, hy 
the method given ahove, ty use of formula (23) an expres- 
sion is arrived at for the stream function \l/ defining 
the flow of gas from a vessel of this kind; there is^ oh- 
tained: 



n 



2^1,= ^6 ~ Ti fjL^^ _£«- sin 2n0 (56) 

All curves \|/ = constant . in the T,9 region start from 
the point T = 0 and meet again at the point T = t^, 
9 = 0. 

The velocity potential hy formula (24) is determined 
hy the relation: 



H n) = C + i f -^-I-, .-(X- t)-P 

2-/ T(l _ T)P 

00 



+ (1 T)-^y i ^ COS 2n9 

^ n y y 

1 ^ n ,0 



With the formula for cp, it is not difficult to 
set up the equation of the jet^ By formulas (8) of part I 
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cos fi _dx .dy sin 0 ^()y, 

whence 

V274f = sin0|-h(l-T)-cose§. 

If use is made of the formulas for <p and ^, there is 
o'btained 

CV3 



~V^2Ta^=— 2(1— t)-' y -^(cos -cose(l-Tr 

1 

Integrating with respect to 9 yields 

u ^ — . . c V/- Y V« rsin(2H+l )& , sin(2n-l)ei ^ 



1 

Since, for 6 = 0, it should follow that y a 0; therefore 
*(T) =s 0, Thus . finally. 



'^L\zJ y,/^\ 2«+l 2«-I J 



1 



1 

The series 



- V/ • \" ^" r ^'"(-»- ^>Q . sin(2»+ l)01 



I 

may be put Into tne form 

rsin C->«-l)0 sm(2»+l)e] 
1^-4 -2^^ZT- + — 2M^1~|' ^^^^ 

where 

By the theorem proved at the end of part II represents 
a series of magnitudes decreasing with n. Thus the expres- 
sion 



Zp, I sin(2n— l)g sm(2n+ l)Q] 
^ I 2n— 1 + 2n^\ J 



and i approaches 1, or, in other words, as T approaches 
Tq, tends to the value 
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VI rsin(2»— 1)0 , an(2»+.l)0 1 
^=L[~2^1 2n+l r 

and the same limit Is approached 1>y tne more oomplleated 
summation (58) entering the formula for the coordinate 7. 
Transforming I into 

CV5 



• h I oV< siP(2n— 1)6 



yields 



05 'X> 



where H indicates that the real part of the expression 
must "be taken. It can "be readily seen that, for a con- 
tinuous change of 9 from 0 to 9, 



1-^cosQ — t'sinQ . % /l-f-cosQ 



and therefore 

CO 

isin(2n— 1)0 



2w— 1 ~ 2' 



1 

depending on the sign of 6; thus 

2 = =4=^ — sin 6. 

In view of the importance of the relation obtained, 
iaore rigorous method of its deviation will 'be preeented. 

Starting from the equation 

m 

8in2m9 



2,cos(2«-l)0 = -2^ 



1 

and integrating it within the limits 0 to 9 results in 

m e 



V< siii(2M— 1)0_ /•siii2iK0 0 . 
L 2»-l J 0 sinO*^"- 

1 0 



tnsiu(2w— 1)0 1,. fsinuO 0 . 

But this limit, as is known, equals ^ if 9 it positive 

and - ^ if it is negative. (See reference 7 > vol* II, 

p* 233*) Hence» 

cc nt 

/=-sm0+22_^— =_an0+21ini^,_^2j— k=T = =^ g 
1 1 
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The second series in equation (57), for T - 
approaches the expression' 

..oo' . 

Y X psi^CSn - 1)9 _ sin (2n + -1)6 
^ ^ 2n ^ 1 2n + 1 

vrhich is a convergent series for all values of 9. OPhis 
can "be shown hy considering the remainder term of this 
series: . 

00 

^ _ V TsinCSn - 1)9 sin(2n + l)9l 
■**n " / ' : — r ~- — I 

— w 2n 1 2n + 1 . 

n 

Substituting in the preceding equation for Xj^ its ex- 
pression given "by equation (49) yields 

E = nr^^2^ 5in(2n - 1) 9 ^ ^^T^nr s^^^^ ^ - ^ ) ^ sin(2n-Hl)9 
^ ^ ^ 2n - .1 .,,?v4i'- ' ZrT- 1 2n + 1 

Inhere k is a finite numher and \^ a proper fraction. 
Hence limj^^Q^I^j^ = 0. 

By now putting t = Tq in the formula for the coor- 
dinate y, the equation of the jet houndary is arrived at 

v^(i-T )Pl!y = :p n^-Yx^ [^injgnj^ - 1)9' 

° C: 2 ^ i- 2n + 1 2n - 1 J 

1 

TT 

where the upper sign of - corresponds to 9 greater 

2 

than zero* 

• «* 

If the width of the infinit ely . distant part of the 
jet is denoted hy -2h, then 



2hy2aTo (1 ~ Tq)P = q 



for ^2a^Q is the velocity at infinity,, -p = Pq(1 ^ t) 



TTy_-^irr 'r'<^sin(2n-l)9 sin(2n+l)9'*!^ fKo\ 

" • i "Z L ""in™ "iin J ^ ^ ^ 

1 

In order to determine the jet contraction, it must 
he noted that, for the conditions under .consideration 

0 < 2 0 V" T^ * contraction occurs at infinity as 
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in the case of the outf Iott of an Incompressible liquid. 

Thus, if the maximum contraction were at a finite dis- 
tance from the orifice, it would then be followed by an 
expansion: The streamlines would be turned by the con- 
cavity toward the inside of the jet; the pressure would 
drop from the surface inward and would reach a minimum 
at a certain point on the line of symmetry. At this 
point the velocity would receive its maximum value, which 
result is impossible. Thus the contraction will be equal 
to the ratio of the width 3b at infinity to the width 
2a of the orifice of the vessel. This ratio is deter- 

TT 

mined from formula (59) by substituting 0 = y —a, 

3 

and taking the upper sign of the first term on the right- 
hand side. Then there is obtained 

00 

TTa _ TT \ 4n 



2b . 2 /j^- 4n2 ^ 1 



whence the contraction is 



b ^ ] Z , (so) 

• CO 

a 



TT + 8/ (--1)''"^^ £ X 

^ 4n2 ^ 1 

1 

CO 

^ ,n— 1 n 



n,o 



The series S = ) (-1) ^ x^ ^ must be. con- 

V - 4n" ^ 1 
vergent since it is an alternating series with numerically 
decreasing terms. Another way of proving it is by substi- 
tuting for Q its expression, formula (49), of the pre- 
ceding section. Thus the remainder term S of the series 
is found in the form 




00 



n 



00 

, n-1 



4 2n-^l ^ xi^^^ 

n . , 

where it is clear that it approaches zero with increasing 
n, for k is a certain constant and \^ a proper frac- 
tion. 
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For Tq = 0, = 1 a^ad this leads to Zirchhoff's 

formula 

a TT + 2 

This will be the approximate value of the contraction for 

OO 

— I f . n— 1 n 

small flo^ velocities. To compute \ k-1) - — ^ tl^ 

L ^ 1 

for finite velocities and finite difference in pressure 
between the reservoir and the medium into which the dis- 
charge occurs, use is made of the approximate formula for 

^n» 

The limits between which the quantity ^ is con- 

a . 

tained are fotxnd "by taking an even or odd number of initial 
terms of the series 



4n2 - 1 



In the first case the positive terms are replaced by the 
lower limiting values, and the negative terms by the upper 

limiting valnes; in the second case the reverse is done. 
This computation will be made for the limiting case 

Tq = — ^ considering only five terns of the series 

3p + 1 

for determining the upper limit and six terms for deter- 
mining the lower. Use is made of the values of , . . . 
computed in part II, and the upper and lower limits 

Zn-*i 
— — ;scn are denoted by A and B, respec— 
. 4n^- I ^ 
tively, to find 

^ _ 0.5254 _ 2^X 0, 47 Q 3 ^ 3J( 0.^4343 

- tJ'^^d^^^ + 5 X 0 , 3 872 _ 6 X 0, 375 5 
63 ^ 99' 143 



whence, with an accuracy of 0,001, 
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B = 0.128 

By rejecting the last term and adding the possible errors, 
A is obtained. The errors will add up only to 0.0003 
and will have no effect on the accuracy desired. Hence 

A = 0. 128 + 0. 015 

Substituting thfese limiting values of the summation in 
formula (60) yields 

TT ^ "b ^ TT 

TT + 1. 14 a TT + 1. 02 

or, if the computation is carried out, 

0.73 < i < 0.75 is approx. 0.74, ! 

a a 

Thus the jjt_ex£anas^ wjj;h__in^^^ 
X.e_S.er5L9jr"Io~ihtrJLtiLLt^ Its' extreme dimensions 

in width are 0.61 X 2a and 0.74 X 2a where 2a as 
before denotesthe width at the orifice. 

An approximate functional formula for the contrac- 
tion is obtained by making use of the approximate expres- 
sion for the function Xj^ given by equation (45), which 
may be changed to 

X, = 1 - pa - !|! E 

. where 

M . ^ 



K = L + 



n + 1 2n(l - Ps) + 3 + O + 2)a 



The coefficients L, M, N do »ot include the parameter n 
and have the values 

L = — ?— " M = - -Jl—l 

1 - Ps 1 + (3p + 2)s 



11= 2(2p + 1) 

(3P. +2)8 



r __^_r_?i! } 1 

L 1 + (3B +2)8 1 - ps -1 
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If, as before, it is assumed that 3 = 2,5, there is 
obtained finally 

^ ~ 2" ~ 4(2— 5s) 2-f 1 9s «-f 1 n-|-u [(2— 5s) - 2+1 9sJ ' 

where 

6-f-!>s 



' 4— lO.s 

The series S entering equation (60) readily can 
be computed by setting 



CO 03 

«2 



«=i(-')"-'4-;^^'-.=['-:r-4iiyi:(-'>"-'4-^^ 

1 1 

(61) 

cv 



2-fiys^^ '(«+l)(4w«— 1)"^ 



cv 



1 

where for simplicity Is replaced "by s. The first 

summation la equation (61) for S Is of the form 

:i .} 5 ' 5 ^ 7 7 J 4 

It is necessary now to return to the commutation of 

.(a)=2;(-i)«-^^3j^ ; 

The particular case of this series corresponding to the 
value p. = 1 will he the second summatloji In the equa^ 
tlon for S, And aiy,) may he expressed in the follow- 
ing form: 

1 

»w=j£(-i)-^:rr'«- 

0 

The series under the integral algn can be summed. Thus 

l^ZT— xp+ 3—3—5+5+7 -7-9 +-J • ■ 

Since the series within the brackets Is absolutely con- 
vergent for t < 1, the order of its terms may be 
changed and this ^ives 

867 
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5 y + ^- l-p- -^arctgV^J 



Thus 



1 3 



4,0.) = 42;(-l)-'j5^:^3j5 = J -j, \l-,).rci&/T<U, 

0 

or, "by integrating by parts. 



(62) 



whence 



4c(l)=| — ^-f-|lg2=0,5437. 



d 3 ■ o 

When the results are added, there Is obtained for the 
contraction 

b ^ 



(63) 



6 ^ 98 

a(|JL) is dttermlned by equation (62) 



where \x = 

4 « lOso 

and 8q by the ratio of the presstire in the Teasel to that 

To 

in the free space, Sq = • which by the formalas of 

1 ^ To 

part I is 

7 

f =(i-Tj—Mi+«o)'-^''-(i+««)' 

Of greatest interest is the Jet contraction for a pres- 
sure near the limiting value - that is, for which the yeloc*- 
ity of the escaping Jet is equal to the Telocity of sound 
propagation in a gas at rest of the same physical state. Thl« 
limiting pressure corresponds, as has been shown, to the Talue 

11 /I 
To = , So = — aai lias the Talue po - 

2|5 + 1 2p 



\ 2P 



- 1.89 p^ ~ that is, 1.89 atmospheres - if the pressure in 
the free medium is equal to atmospheric. If 

^ 4-105, 

ie set in formula » Jl. 0«2, the result is fx » 3. 9. 

To compute accurately the definite integral, 
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J 1 + 



for such value of i-L is rather la'borious. In view of the 
f act , towever , that this entire computation is of an approx- 
imate character, the problem may "be simplified: namelj^, 
compute j(4) and j(3*75) and then hecause of the near- 
ness of the values of these integrals, find j(3,9) by 
simple interpolation (assuming proportionality "between the 
increment of the function and of the independent variable). 
There is readily obtained 

-1 • • 

J (4) = / «i dt = l--i+i^lg3 = 0,1402 

J 1 t 2 3 

0 

1 11/4 _ 

j(3.75) = r 1 dt = 4 (1 ^ i +L\lL+y2 Ig cotg H- 

v/ 1 + t \3 7 11/ 72 ® 

rz 0.1500 

whence j(3.9) = 0.144; and by equation (62) 

4a(3.,9) = -.2 TT_ + 15^ — 0.144 =.0.212* 

8. 8 8. 8 X 6.8 8.8 X 6.8 

Substituting this value of a( jJi) and the correspond- 
ing So in formula (63) yields 

8S =^ 0.5 + ^ X 1. 087 + X 0.424 = 1,08 (64) 

29 145 

the mean value between the limiting values of the series 
8S obtained. The contraction then is given by 



a TT + 1. 08 



♦The procedure for checking is as follows: Compute 

the accurate values 4a(4) = 0,2079 and 4a' ( 3. 75 ) = 0, 2193 

whence by interpolation again it is seen that 4ct{3.9) = 

0.2124, a value agreeing with that already obtained. 
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This coefficient decreases with decrease In the pres- 
sure In the reservoir "because of the increase in value of S, 
This decrease is sufficiently uniforn«, as may "be seen from 
the table: 





0,2 


0,182 


0,154 


0,137 


0,117 


p, 


1,89 


1,79 


1,65 


1,5U 


1,48 


h 

a 


0,74 


0,73 


0,71 


0,70 


0,G8 




3,90 


3,50 


3,00 


2,75 


2,50 



(65) 



Finally, the expression for the quantity of outflowing 
gas is given by 

£=2a.^V^7?,(l-hT,)^ 
vhere is the density of the gas in the vessel (at a 

far-removed region from the orifice) 2a, as before, is the 
vidth of the orifice and a is defined by the formula 



where 



Y=l,40=l+A. 



It is necessary first to consider the case of the 
outflow of a gas from reservoirs with various pressures 
into a space where the p re ssure _i^ con stant (e.g., into 
the atmosphere). Then "* 



L 



1 



and the final formula 
E 



y = 1.40, P . 2,5 



OP. again, If V = 1.40, p « 2.5 



(66) 



(66') 



HACA TM ITo. 1063 



63 



Since the contraction - is a function of and 

a p o 

therefore depends on the ratio --^ , V7ith constancy of 

the ratio, the discharge quantity of the gas is propor- 
tional to the square root of the density or inversely 
proportional to the square root of the temperature. 

THovr suppose that the state of the gas in the reser- 
voir remains unchanged and consider the flow into a me- 
dium of varying pressure. The velocity of sound in the 
gas corresponding to the same physical conditions at a 
great distance from the orifice is denoted 'oy Cq, Then 

and the formula for the discharge may "be given as 



or, sTibst itut ing y = 1^4, 



/ 10" " ■' ' " 2 ■ 

In this formula for E only the last two factors 
that depend entirely on the pressure ratio vary. 

As regards the jet contraction or the dijohargo 
coeffici ent as the magnitude h/a also will be called,, 
it is accurately determined by formula (sO) and approx- 
imately by formulas (63). For an approximation of accu- 
racy up to 0.01, numerical values have been given of 
this coefficient for the limiting pressure and several 
other smaller pressureso It is found that a practical 
result of the s?ime accuracy is obtained if, instead of 
the cumter some . formula (63), the following expression 

b 

is taken for =: X; 

a 
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•Vi ■ TT 

X = - = (68) 

T7here . k .is a certain constant. This formula gives very 
good values of tho c ontract ion. f or small difference in 
the pressures 'between reservoir and free medium - that is, 
for s near zero; for Sq = 0 the values of b/ a "by 
(60) and (68) a^ree \7ith the Eirchhoff formula for the 
Jet contraction. The constant k is so determined that 

for s = — = 0.2 the results "by (6S) and (53) agree* 
?or this it is necessary that 

(2- ~ k)(0, 2 - 8S(C.2)) 

and, since .83 has been found (in, equation (34)) equal 
to 1.08, k rt 4.6. Tho values of 8S for the values 
of Sq given in table (55) are correspondingly equal to 
1.08; 1.15; ] .28; 1,35; 1.44, To these correspond 

1.03; 1,17; 1.2&; 1.37;. 1.45 

numerical values of tho binomial 2 - 4.5so entering 
formula (68) in place of 8S. 

Such difference has no effect on the results for the 
second decimal accuracy Trhich has been assumed. 

The agreement v/ill bo even better if the discharge 
coefficient is expressed by the formula 

% ^^0 ^ !: ^ (69) 

^ TT + 2 - 5s + 2So 

The series of values of the function 2 - 5sp 2so 
the same Sq vrill, to an accuracy of 0.01., be equal 
1.08; I.IS; 1.28; 1.35; ?s.,nd 1*44, rrhich are equal, 
respectively, to the above obtained approximate values 
of 'as. entering the exact c ontract ion formula. Thus, 
formula (",2) or (6S) for k - 4.6 quite uell expresses 
the function \ and may very conveniently be applied for 
practical purposes. The ,di scharge formula (57) for the 
assumed round values of y - 1,4 and g = l/(7 l) = 2.5 
tak a sufficiently simple form if the variable Sq is 
introduced in it. The latter, as has been shoirn, is con- 
nected with the pressure ratio Pq/Pi* 
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(70) 



By transforming equation (67) in this manner, the discharge 
formula is obtained as followss 

'^=2ac,,p;A^3^°^^^ (71) 

It is noted that of the two simplified formulas for \ 
the second, (69) , is the more rational. Thus, to turn for 
a moment to formula (60) 

x = 1 

and substitute in it for ^ the trinomial 

expresses approximately the function Xj^. From the equa- 
tion for Xjj 

which gives for Xjj( 0) = 1, xi^(o) = -^ 

Suljstitute these values in the above trinomial and compute 
the series S, entering the denominator of the formula for 

but 

1 I 

" 0 
or carry out the integration 

1 

(4>7.' l)'(H-f-l) 3 



whence 

Assume as before p = 2.5, which givej? 



and 
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Limiting to an accuracy of 0* 01, cons ider this formula 
identical with (69). 

It is considered of interest, finally, to call the 
reader's attention to a very simple connection between 
the variable Sq by v/hich all the characteristic con- 
stants of the problem and the temperature of the jet is 
expressed. For the density and the pressure within the 
gas flow the formulas are: 

hence 

P Po Po 



and since by Mariotte and Gay-Lussac's law p/p = RT, 
where T is the absolute temperature at the point con- 
sidered, the foregoing equation may be rewritten a^ 




since s = — ^ — (see pt, II), Applying this relation to 
1-T 

the part of the jet remote from the orifice and denoting 
the teraperature of the gas there by yields 



"^0 = 



3^0- Ti 



(710 



This investigation on the outflow of gases will be 
supplemented by comparing the results obtained on the one 
hand with approximate theoretical formulas applied for 
computing the discharge and on the other hand with the 
results of tests, Puroly empirical formulas are not dealt 
with although some of the latter well express the phenomenon 
within certain limits, as, for example, the formula of 
Parenty (reference 8), 
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^'or a rational basis of the approximate theoretical 

treatment the adiatatic lav/ was assumed (also in this 
investigation) the assumption being made that the out- 
flowing jet at a certain distance from the orifice has 
the maximum contraction and that at the points of this 
contraction the velocity of the gas particles is constant, 
As a result the following formula is obtained for the 
discharge formula: 




where S is the orifice area and X the jill^harge 

c oef f icj>ent equal to the ratio of the area of the con- 
tracted cross-section to the area of the orifice. 

The above equation does not differ in form from 
equation (67). the only difference being that the dis- 
charge coefficient was not determined for any, or even 
for a particular shape of orifice. It has usually been 
assumed that it has a constant value depending only on 
the. shape of the vessel and orifice. Such as suu^pt ion, 
as is seen from the problem solved here, is far from 
true. In this case this coefficient, for a change in 

--^ from 1 to the limiting value 0.53, increases from 

OoBl to 0.74c The increment thus cor.stitutes more than 
21 percent of the lower limiting value. If the orifice 
were round and not in the form of a slit, as in this case, 
a still sharper difference in the values of \ should be 
expected, for then the lines of flow would converge toward 
the orifice from all azimuths and not from two as is true 
in the present case. For this reason, wheh ,it was attempted 
to 'apply the discharge formula v/ith constant K to the de- 
termination of the true discharge, experiment did not turn 
out to be in agreement v/ith the theory^ In view of this 
'Parbnty (reference 8) relying on the tests of Hirn (refer- 
ence 9) assumed that to apply the formulas based on the 
adiabatic lavr of pressure " change t^o flow ' discharges from 
orifices v/as incorrect. Hov/ever the results of Hirn»s 
tests v/hich ho presents show precisely the increase in the 
discharge coefficient X which is predicted by the present 
theory^ The possibility of such a variation v/as foreseen 
by Parenty but having remarked on it gave It no further con- 
sideration since he had no means of making a quantitative 
estimate of the increase in X. 
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Another fact is considered here that is of interest. 
Having obtained the discharge formula (67a) Saint-Venant 
called attention to the following paradox. If this for- 
mula vrere applied for any ratio of pressures in the res-- 
ervoir and the open medium, the discharge, increasing from 
zero, would pass through a maximum for a certain pressure 
ratio and thereafter should again deer eas e, "becoming zero 

Pi/Po 0' value of P for p^/Po corresponding 

to the maximum discharge is determined from the condition 



1 y-^i 1 



hence 




'This is just the limiting pressure ratio correspond- 
ing to the instant when the gas in the contracted part 
of the jet moves v/ith the velocity of sound propagation 
at that point, as remarked hy Hougcniot (reference- 10) . 
This condit ion /Cannot of course occur in practice. When 
S.?int-Yenant conducted his tests on the flov; of gases he 
found that on lov/ering the pressure in the free medium 
and varying the ratio, Vi/Vq ^rom 1 to 0.53 the discharge 
increases; hut on further lowering p^ the jjrocess becomes 
regular, there being no further increase in the discharge. 
This surprising result v/as lon.g looked upon v/ith doubt but 
Hirn's tests, conducted not very long ago, confirmed the 
results of Saint— Venant with the difference, however, that 
Hirn observed an increase- in the discharge beyond the limit 
indicated by his predecessor. According to Hirn»s tests 
the discharge reaches the maximum value for P^/Po = 0.26, 
approximately. The change in discharge on lowering Pi/Pq 
from 0.53 to 0.26 is, hov/ever, ins ignif icant, f or which 
reason this may not have been noted by Saint Venant in his 
less detailed observations. 

If it. was attempted to apply formula (70) for \ 
determined by relation (69) beyond the proper limits of 
its applicability the same paradoxical result wotild be 
obtained except that the maximum discharge would corre-*^- 
spend to a value of Pi/Po somewhat less than 0.27^ a 
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value very close to Hirn's limit. This interesting agree-*- 

ment shovs that the' present formula expresses sufficiently 
well the investigated phenomenon in its essential features. 

It is now natural to inquire into the character of 
the motion in the case where the pressure in the free 
medium into which the jet discharges is lower than the 

limiting; that is, if ™ < 0,53, If it is assuxaed that 

Po 

flcvr x-emains steady with continuous change in the velocity 
and pressure v,rithin the "boundaries of the movin^R: gas uiass 
the region of the variables g would De in the form of 

a seiulcirole of radius > ^-J- — ^ This is the very region. 

2 3 1 

considered in part I where it v;as shov/n that steady motion 
of the type that is of interest to us was not included in 
the number of possible motions, Hougoniot ( reference 10), 
• states the followin^:? i- p^. < 0,53 p^ the escaping jet 
is divided "by the surface' over which the velocity of the 
particles is equal to the velocity of sound, into two parts, 
the pressure in passing through this surface changing dis— 
continuously; c^.bove this partition surface in the Jet the 
pressure is equal to 0,53 p^ and below it is equal to p^^ 
(This phenomenon reminds Parenty of the separation from 
solid bodies,) But the flow of the gas is considered as 
steady in both parts of the escaping jet and the surface 
of pressure discontinuity as everywhere normal to the 
streamlines. 

This latter supposition appears highly improbable 
since the character of the motion in the upper part of 
the flow should radically change immediately after the 
pressure in the free medium passes beyond .the limit of 
0.53 Pq. In fact, first, of all it can be easily shown 
that the width of each elementary tube of flow will be 
a minimum at the point where the limiting pressure occiirs. 
This is because the cross-section is determined as the 
ratio of the quantity of gas carried by the tube divided 
r- P 

V So* ( '^Q ) and this denominator passes through 

a maximum at T = —i---, • Therefore taking the tubes of flow 

2g+l ^ 

normal to the line T = -2^-^ a Eiinimum discharge of ^^as 

from the vessel shall be obtained for the case where this 
line is a segment enclosing the orifice. The discharge 
coefficient will then evidently equal 1 and, therefore ' 
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in passing through the limiting pressure in the receiver 

this coefficient and the gas discharge should immediately 
increase "by more than 30 percent, a condition that is in 
entire disagreement with reality. 

It is assumed on the other hand that the phtno^enon 
could ho explained in the following way. Together with 
the authors referred to it is- supposed that the jet is 
divided "by a certain boundary surface on passing through 
which the pressure changes very sharply. It may he 
imagined that the trace of this surface on a plane parallel 
to the flow as a curve supported at the edge of the orifice 
further on the curve resembling the ciontour of a tongue of 
flame moves into the open medium, Ahove this limit (inside 
the vessel and. the adjoining part of the jet) the flow will 
he stable and the pressure drops from p at the far re- 
moved parts of the vessel to 0,53p^ on the described 
boundary curve,- At the remote part, however, the jet 
forms waves,* Tiiese waves have an enveloping boundary 
curve. In a very thin layer of this part of the jet ad— 
joining the curve the mean pressure will be < 0.53Pq and . 
the velocity of propagation of sound c^ < being the 
same velocity for the boundary layer lying beyond the 
botindary curve. The lov/ered pressure tends to be prop— 
agated beyond the boundary curve, following along the jet 
in the form of a plane wave. But this v^ave is carried 
backv/ard by each infinitely thin jet element and since 
the velocity of the gas particles is also c^ no waves 
are observed in the upper parts of the gas flow,. In order 
that the boundary curves may serve as an envelope of the 
waves approaching it, it is sufficient, as it appears to 
assume that the velocity of the waves normal to this line 
is the same v;hether the wave moves upward or downward. If 
\ denotes the angle formed by a jet element with tha bound- 
ary carve passing through a given point, then having deter- 
mined both normal velocities by Hiemann's r\xle equating them 
and applying very simple hydrodynamic considerations there 
is obtained ^ ^ ^ ^ 

*These waves have been observed and studied recently 
by Emden, The results of his tests are described 'in refer- 
ence 4» The waves appear immediately ' aft er the pressure 
in the reservoir drops below 0,53pq; their length increases 
v/ith the lowering of the pressure in the reservoir. Smden 
also gives a theory of the phenomenon v/hich, however, is 
entirely unfounded. It is sufficient to say that notwith--* 
staading the existence of waves Emden considers the pressure 
throughout the jet as constant, which of course is impossible, 
and makes this assumption the basis of his analysis. 
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sin\ = 




If the pressure and density varied, even very sharply^ 
"but c ont inuously , p^^ ; =^ 1 , Thus it is seen that 

^t. £.i£i£ iSl^XlSct^ iiiS i^i^^SrZZ si a con-* 

V/ith this hypothesis steady flow above the houndary 
curve may be determined strictly luatheinat ical regardless 
of v/hat occurs in the remaining part of the jet^ It is 
not difficult to show; namely^j that pn this curve, given 

by equation t = — ^ — - the relation holds 

^ 2B+1 



cptgX(l-.Tj^-± = Q 

one sign corresponding to the left half of the boundary 
curve, the other to the right half and x]/ and cp denoting, 
as before,, the stream function 4nd velocity iDotential. In— 
eluding this relation among the boundary conditions it may 
be shown next that, together with the other conditions, it 
is entirely sufficient for the determination of co and 
in the 9 region. Having found cp and ^ it is easy 

to determine the gas discharge per second^ It appears that 
if this discharge v/ere strictly constant or changing slightly 
with change in \ from zero to its limiting value, the 
explanation just given would be near the truth. Incidently 
it may be said that the limits within which X may vary are 
not v/ide; this angle will not be large^ Por this reason 
the relation previously given between ^^ and cp in all 
probability will give a result not deviating too much Irom 
that which would be obtained by simply taking 
along the boundary cur^'-e^ A small variation in the dis^ 
cha.rg8 may also be expected from the consideration that 
its value will depend on /sin\ds extended over the 

boundary curve. This integral is evidently equal to the 
total length of this curve multiplied by sinX, and its 
length will decrease with increasing \. It may be noted, 
finally, that in assuming the above explanation of the flow 
phenoiiiena there is obtained an entirely continuous transition 
from the problem solved above to those cases v/here the given 
analysis is inapplicable^, A mathematical treatment of the 
proposed hypothesis is intended in the near future. 
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PART IV 



PHESSURE OF A GAS JET OIJ A PLATE 



The study of pressure of a gas jet on a plate will- 
tegin T/ith the consideration of the impact of a gas jet 
on Si plate perpendicular to the initial direction of the 
jet, assuming that the jet is symmetrically divided into 
two parts hy the plate* Again the corresponding problem 
for the case of an incompresgihle liquid is used. The 
solution of this prohleiu is given in the paper hy 
JoukOT7sky (reference 2), By use of the same variables 
as in part III , 



^1 + i'^i 



z = X + iy , ig Vq = i& + ie 



dw 



sin^O - U^) 



(71) 



where m is the angle of inclination with the X axis at 
distant points of the two parts into which the jet is 
divided by the plate. For the X axis the line of symme- 
try of the jet is taken, the Initial direction of the 
jet being in this case parallel to the X axis. 

It will be sho-vn that formula (71) expresses pre- 
cisely the required liquid flo^?, ■ Attention will be 
directed first to the range of complex variables w and 
^ + iG which correspond to the flow sketched in figure 
5, The region w is "bounded by two str^^iight lines par- 
allel to the real axis, symmetrically placed with respect 

■to it at a distance In the sketch the outer bound- 

2 • • 

arios of the jet EA^ and DC* correspond to these straight 
lines. The flow boundaries C30 and ,AYO correspond to the 
upper and lower sides- of the positive part of the real • 
axis of the region w; ^the point 0 corresponds to w = 0. 



The region f + iQ . is bounded, in the first place, 
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"by the segments of the straight lines parallel %o the 

TT TT 

axis and having coordinates 9 == o and 9 = - - lying 

*to the right of the imaginary azis;, in the second place, 
hy the segment of the imaginary^ axi s" lying hetwe^^^i the 
ahoverinent ioned parallels. On the sketch the straight 

• Tf 

line 9 = ~ corresponds; t o the right-hand part - of the., 
plane and the line e = - 2 to the left-rhand part. The 
segment of the imaginary axis included "between the points 

TT ' • 

9 -='^ and 9 = m corresponds to the "boundary YA of the 

flow, for here the velocity is v^,^ = 0. The segment 
syipmetr-ica,X • 1 0 that, just mentioned corresponds to the 
curve BC« Finally, the "botindary CD is represented in 
the d ie region by the* segment, included between the 
points; 9 = - m and 9 := 0 . of the imaginary axis and the 
curve EA* by the segment bounded^ by the points Q ^ 0 
and. 9 = m, ■ . . ^ 

l«'ow , proceeding along the. boundaries of the . -9 + i9 
region, the author will show that the point w will then 
describe the above-mentioned boundary of the w region, 

With the point 9 =: as the starting place, it can 

be seen, from equation (71) that for these values of d 
and 9 w = 2kiTi , where k is an arbitrary integer - it 
will be taken equal to zero. If now the point -6 + i9 

moves along .the line 9 ^= then w moves along its 

real axis at the upper side of this axis, since for 9 = 

^ - €, then, at infinitely small 'J^i = k €, where k 

is some positive quantity. When the point -6 + i9 arrives 

TT 

at the position ^ = 0, 9 = w will have -passed along 

the segment of the real axis from 0 to 9i =? ^ lg(.l - sin^m) , 
As + i9 moves farther along the imaginary axis, the 
point w will continue its motion along the rp^ axis in. * 
the same direction up to cp^^ = <» corresponding to -6 = 0, 
9 = m. In passing through the point = 0, 9 = m-, the 
logarithm in formula (71) receives an increment - Tri and 
for iS=;0, m>9>0, w will move forward along the 
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straight line w = i~ from +<»Hr9' tcH^<»t-i2.; this position 

:7 corresponds to 9 = 0. In passing through the point 
^ = 0, e = 0, the logarithm of fornnila (71) \7ill receive 
an increment Sni and vr will pass discontinuously from 
the upper boundary of its region to the lower and. will 

^ ^ 0 0 . « 

move along it from to + 00^ i^ as ^ ^ xB 

^ 2 

moves from 0 to - i m, Further, as + i6 moves 
through -im, the logarithm increases hy —Tri ^ w jumps 
to the point + 00 and along the lower side of the posi- 
tive ^art of the axis returns to its initial posi- 

TT 

tion as ^ + i9 moves from -im to ^ 2 from 

— i 21 to 00 ^ i U along the "boundaries of its region. 
3 2 

Thus the fact is shown that formula (71) is an ac- 
tual solution of the problem of the impact of a liquid 
stream on a plate. In order to solve this problem for 
the gas jet, it is necessary to proceed according to the 
rule given. The expression w is expanded into .a series 
and its imaginary part separated. 

' Thus , 

2 Ig sin tJL 

Q ^ ■ i 

^= 2 Ig sin LLJl ^ Ig (cos2m - cossl-t-M^ lg2 

. ^ + i9 T . 1^ + i9 -im . ^ ^ iS i9 + i m 
= 2 Ig sm Jg sm Ig sin : 

i i i 

Introducing exponential functions in place of the trigo- 
nometrical transforms this into 



^ Ig (^^iri^ tJLJl ^ sin^ m) 
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i-osin(0-|-iiM— isin2»(0^->»)— 2i-os2n04--2/si 
or, after reducing. 



= — - 2,. — J7~ (1— co2wiw)(cos2«0— /siir2wO), 



1 

end, finally, for \|f ^ 



= — - — sii>2>/;(l — ros2;^v>n. 



I 



Since (v =yF^) 



0 



^ = 2^ ^ ^ ~ ^ sin2«0(l— cos2nm). 



This is the expression of the stream function in 
the variahles T , 6 for the liquid flow^ Hence for the 
gas flow the same prohlem should be solred by the formula 



f=Z7,{f.)"j7/°2««(i-c«»2««,. ^^,) 

1 

According to formula (24) of part I the following 
expression is obtained for cp: 



111 = C-(l-t)-^ y^-^ x,cos2«0(l-cos2«,«). (73) 
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If it is assumed that C = 0, then, as in the case 
of the liquid flo^jr, T= 9 = 0, cp = 0. and \J/ = 0, 

The series for cp and \|/ are ahsclutely and uni- 
formly convergent for any T < Tq as is cle^.r from the 
general considerations of part II, and therefore the 
problem is solved hy the relations (72) and (73), 

The vridth of the plate and the pressure on it ^ill 
nc:7 be determined. The width of the plate uill be de-- 
noted by 2 I, Then I is found by substituting in the 

TT 

expression for the coordinate y the values 9 = - , 



But 



by by 
dy = — dcf + — d\|/ , 

Sep S\l/ . 



6-2- 



and, since for 6 = H , \|/ = 0 = constant 

2 

8 = n 

(dy) = ^ dp = ^ dT = / ^ |S dT 

Sep bcp dT / J^ar 

6=3 

By su^bstituting the formula oTjtained for cp and noting 
that 

I = 

there is obtained 

CD '''o 

TT /2a. nT 0^ yn , 0 J dT 

.n+ 1 



1^1 yO (74) 



In conputing the integrals entering this series, it is 
noted first of all that 
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and on the basis of the equation for Zj^ 
Hence 

0 

To 

0 



On the other hand 



0 



0 

and BXisce on the 1)a8l8 of the differential equation for 
the function z,^ 

(1— '»"~n- rft ' »" 

therefore 

0 ^ 0 

Comparing this relation with (75) results in 



or 



,Since = Tny^. ^ + 5 T;'' " =^n. 



1 

a 



2m . . - 
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With the aid of formula (76 • ) , equation (74) can "be 



written in the form 



(77) 



1 

ihe first of the summations entering the a'bove formula 

is an absolutely convergent series. It is computed as 
f ollovs : 



1 1 

CO 



-|h-J+..-)=|(§-i)- 

cos2//wi ri . cos(2w— l)iw 



1 

It is 7e«dil7 ■••n« liow«Ter, that 



tr V / 1 c*><(2tt — l)w» — *siD('2ii — l)»i 



1 ess 



=J 2] (— i)""' lcos(2M-ni»»-isin(2»»— !)»»>// = 



0 1 

1 <x 



0 I 



HAOA TU ITo. 1063 



Since Id ^H^i ~ * Integral obtained 

1 

is expressed as follows: 

1 1 



(areegfe— * = i arctg.-"- j J ^ii^ 

t\ A I 



0 

1 



^arctge __ + _J__ = 

0 ' 

1 

= "2 (l+r'"')arctge— - ^ ; 



and since 



therefore 



r. / , i /14-sinw? 

^ (1+C0S2,,,)- + Ig 1+^5??^ + 

^ \ 8 ^ 2 ^ 8 1— sinm / 

Heturoiiiig to the initial formula for IL , it is seen 
that the first row of its nev expression gives the expres- 



sion In finite form of the function V (— (^oH^n l)m 

mm^ 4n' — 1 

1 

and the coefficient of 1 in the second row is equal to 



_2;(-ir-'^?^^gEr-; then (78) yields 



1 



1 

and, finally. 

CO 



2:(-i)"-5^=|(i«— i) 



(-1)" >-^~_ = -^l-cosiwj. (79) 

Making use of this formula, equation (77) is finally 
transformed into the following relation 
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(80) 



+2 (-I)""' 1^ (l-cos2»m)}(l-To)-<^. 



1 

The series In the above formula must "be convergent. In 
order to prove this the remainder term Is set up, making 
use of formula (49) of part II for x^. There results 



i?^=»l=:Wo£ (-1)""' |4^-i (l-cos2nm) 4- 



It is clear that th« li«a:aoD " ^^^^^ ^n' 

are proper fractions* 

By passing to the computation of R — the resultant 
pressure on the plate - the pressure behind it is denoted 
hy Pi and it is noted that 

ii=2 J(p-p.)dy=2j {-i i-i^}-2p,h 

0 0 

But 

j^*p'==fc?;'(l-T)^+', 8iJice p=p,(l-T)^ 

hence the integral entering the expression for H — it 
is denoted "by T - on substituting the value of the 
function (p, becomes 



CO 

7^ yj 1-C0S2WM/ 



0 

Further, using the expression 

T 

^W-j|,(l-T)-^(x"-/. + "4-V.)^|, 
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there Is o'btained 



I) 



B7 (76) 

Su'bstitutlng this expression in the sum to lie eo^sputed, 
7i«ld« 

«/2i ^«^.o ^'^^ 



hence 



W2a 

1 

The first of these terms on the hasis of (74) gives the 
■agnltude 

and therefore making use of (79) and noting that 

*?o''(i+^)=;;ifT?o''=«?«. 

t 

there Is o'btained 

r=p,i + ^?-3o(i_cosm). 
The a'oove formula for E yields 

R=Q?j2^,(\-cQsm). (81) 
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By making use of equation (80) and considering the relation 
*/ 2aT . = - the velocity of the jet at very distant 

points, po(l " '^o) " pi density of gas at the 

same points, the required formula for E is obtained: 

E = 2lp, v/ ~ (82) 

+ 8 y 1 ^nxn^o^ (i „ cos ^2nm) 

1 — cosm^ 4n — 1 

The angle m may "be determined "by equation (80) in 
which all magnitudes except m are given. Thus 



Y - 1 y - 1 V - 1 1 -T^ 



where p^ and are the pressure and density at the 

critical point of the "branching line of flow; whence 



1 



Finally, the difference "betvjeen the values of \i/ on the 
jet "boimclaries (a' magnitude will "be denoted "by Q) is 
det err/iined from the condition 

qpo = 2pihVo = 2"bVo(l - Tq) p^^ 

where 2h is the width of the jet at infinity. The mag- 
nitudes h , Vq, and Pi should, of course, "be consid- 
ered as given. 

The resultant pressure, after m has been found, is, 
of course, most simply computed "by formula (81). It may 
he renarlzed here, incidentally, that this formula may he 
derived very simply from the momentum theorem. It v/ill 
then "be found that the formula in no way depends either on 
the shape of the plate, pr ovided , that the latter is sym- 
metrical v/ith respect to the center line of flow, or on 
the relation between the pressure and the density. Thus, 
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denote M the momentum of the gas enclosed in the jet 

"bounded on one side hy the perpendicular section of the 
initial jet at a very large distance and on the other side 
by tv/o siiuilar sections passing through the distant points 
of the branched jet. The increment in M in the infini-^ 
tesimal tiae interval Lt is evidently equal to 
2TDpi VQAt(cos m - 1)vq = Qrpo^o^^^® ^ " since 215 

is the Initial width of the jet and ftp ^ is the quantity 

of gas passing through its cross section per second. The 
impulse of the external forces is given by -EAt. By the 
momentum theorem there is obtained 



and this eqxiation, after simplification, leads to formula 
(81). For this parpose the less general formula (82) is 
of more inportance. By using it the second fundamental 
probien of the investigation may he solved agpr ox i^ma t^e 1^ ; 
namely, the pressure of a boundless gaseous fluid on an 
obstructing plate* 

The approximate summation of the series is started 
by entering the denominator of formula (82). Por this, 
the approximate expressions for the functions x^^ 

already used may be used again in the problem of the out- 
flow of a gas from a vessel (in deriving formula (61): 



ftp ^Vp (cos m l)At = -EAt 



2 



£5sf__ ^ 3 0sf 1 

3(? " 5s) 2 + 19s n -I* 1 




where 



6 + 9s 



s = 



T 



4 ^ lOs 



T 



Denoting for briefness the computed Z by L yields 
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03 

The series of the first, second, and third row will "be 
denoted Hi, Hg, H3 , respectively, and it will "be 
noted that after computing IS^ and H^, is ol^talned 

"by suhstituting m» = 1 i^^ H3. Then can he con- 

sidered as the limiting value of the magnitude Vx'* 

There is ohtained: 



Z/ 1 . 1 \ \^ cos2nm 
(-i)"-(s=r+iis+i)-Z'-'> 



11 = 1 



The first of these sums has the value 
tha sennnd maj be glren In th« form 

I. A- 



(-1)"- --^l^JIT '"""1 ~ 2n-l 



1 1 
and the third in the fora 

Z„ ,cos(2«~l)w' V, ^. 

(-1) -IT-^T" ~ «"'"L(-^) 



. , sin(2n— l)iM 
2n— 1 



2 " 
By the use of these formulas there is obtained for 49^* 
the expression 

/ 
1 
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whence 

^» = "r k=^L^~^^ 2n-l ' 

1 

2^ (_i)-.!i!ii|z-i]i'=|' 



But 



1 <• 1 

* A- 



hence 



2(_l)'«-'cos(2n-l)«.= ^3^J](-l)"-(cos-2««,-l- 
1 1 

+ eos(2n-2),»)= J- [l+(-l)»-Vos2A«» J, 



0 

The last term la zero, for 

tn 



Jeos2hn C cos(2A: — l)wcosm~ siD(2A;— 1 )>nsiD»i 



0 0 



sm(2A:— l)m /' msin m siD(2A:— l)w 

2*— 1 J cosw ' 

o 

and from the properties ot the 7ourier integrals 

m 

C msinm sin(2A — l)m , 
J cosm w 



Hence 



Thus 



cosw* 

0 



-V,=X (-1)"-' jj^-f (l-«.s-.«„.) = -^- I — = 

1 0 

(84) 



"^See reference 7, p, 233. 
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Hov tvxn to the computation of H 3. 



(n-f a)(4n-— 1) 



1 

1 ^i) 



1 (-1)""' ^?=T (l-«>«2«'»)'«. 



0 1 

The function under the integral may he expressed in 
finite form through the lower transcendentals • ITor this 
purpose > consider it as the real part of 

2 (-')-' K ('-«*"")• 

1 

The integral may then be expressed as 

1 GO 

0 1 

where & is u««d to indicate the real part of the com- 
plex q.tiantltjr folloving it. But 



Thus 



1 3 



4]S\=Bjt ^ { (e~*'"—te"")ardg //e ""—(1 — t)(nctg^i ] dt. 

0 

appl/iag integration hy parts, there is o1>talned 



1 r. 



rr / n 
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g— wi gM» 4C0SW 8u.tsinw» 



1 ' 1 io^l 4a''— 1 ' 



thus 



wj -iuisinmr .... (r. mW 

1 1 



0 0 

The Integral in the formula for N3 is equal to 5^ for 
m » 0. Finally, there la obtained 

r(cosiw— 1) , 4uLsium m\, 4{a H'^-^di 



1 



4(x f (l-f-<cos 2»t)f' 



1-f 2teos2m4-i' 



dt; 



4 V t: . 4a sinm , . /- »A 

(85) 



1 

4a-— I ~.) ( 1 -f <Xl +«-'-|r2<cos2«t) 



8a(l-} -cosm) r__JJ\l--Ottt 



0 

Setting in this formula |a • 1 yields finally the ralue 
4Va (Vs is the second suamation entering the formula (83)). 

4.V., T. , 4 s'mm . . /~ ni\ 

T=^n = - 3 + 3 'S^'^U - 2 

(86) 

1 

8(l+cosw) r t{l—t)dt . 

3 J (l+0(l+^'+2teos2?«) ' 

0 

The definite integral in the ahore equation is easily coift- 
putedt it is equal to 

mcoigtu -\- Igcosin— 1^2. (87) 
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Everything required for the computation of 2L by 
formula (83) has been developed and the expression for the 
force on a plate for any gas jet can be set up. Now, con- 
sider the cas6 where the jet is infinitely wide - that is, 
the problem of the action of a boundless gas stream on a 
plate* For this condition m - 0, since the flow after 
passing* roimd the plate must finally resume its Initial 
direction* 

Therefore compute 

lim * ; 8L^ 

m = 0 1 — cosm 

where the exact value of L is given by the series enter- 
ing the denominator of formula (32) aud the approximate 
value by relation (83). Using formulas (84), (85), (86), 
and (87), there is obtained 

lim tl = t, = ^Ig^ r^^^ =; 0.32907 

m 0 1 - cosm 2* m « 0 1 - cosm 3 12 

m = 0 "™coIm ^ 4ijl2 - 1 ^ Tl4\x^ Z'^l) 4|xS Z ij U + 
and s ince 

/ (1 + t)^ 4 / i + t 

'-'0 •-'o 

therefore 

lim = ^-II^* ^ 

m = 0 1 - cosm 4^.^ - 1/ 1 -k t • 2ii + 1 4(4n® - 1) 

0 

Hext, T)j' (83) 
lim = 4 - 10s - ^ --^^ 1.3163 



m » 0 1 — cosm 2 - 5s 2 + 19s 

3 1 1 



+ 



L (2 - 5s) 2 + 19sJiMm, - 1 / 1 + 

iii ^ = P(s) 

2u + 1 4U.S - 1 / 



dt 
t 



HA OA TU ITo. 1063 



89 



where 



I, = ^JLll, s = . s = So (88») 

4 - 10s 1 - T 



Formula (32) takes the forxn 



R = SipiVo^ (88) 
rt + V 



Sq heing deterained by fornnila (82'): 
which, nay be rewritten as 



s 



0 



5c 



3 



"by svilDs t itut ing the velocity of sound c at the distant 
points of the flow and the value 1,40 for Y . 

Por s = 0 the expression found for R gives the 
formula of Kirchhoff (in reference 2); 

S = 2lpaVo3 — - — = 0,44 X 2lpiVo^ 

7T + 4 

applicable for the pressure of the flow of an incompress- 
ible liquid. This value is approached by the accurate 
formila (32) for s = 0, The greater the velocity, how- 

ever, the greater the magnitude — — and therefore the 

TT + P 

reaction of the gas flow with increasing velocity increases 
somev/hat more energetically than in the case of liquid 

flow. How, compute the coefficient — i or the limit- 

TT + P 

ing value of s equal to 0,2 and for values near the 
limit, to an accuracy of 0^01, 
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1 

Setting fl^ » 0.2 yields « 3.9, J ^^^<«=;(jx)=0,U4; 

i>(0,2)=l+ ^ 0,857 + ^ . 0,«9=2,28. 

i: 



, =0,58; 



- the Telocity of sound at the distant points of the 

gas jet; if Pi, the pressure in t:iat region is eq.ual to 
1 atmosphere; then in the ease of air e » 333 meters per 
second* 

For 0,1818; 0=3,5; i;^=cy/jj= ah out 318%,.; 

^u)=0,162.5; &(pL)=^0,508; 
W1818)=4~l,818- 2^ + 3^^ ^ ^ 

3 1 818" ^ 

0,102=2,42; ^=0,56. 

'or ^ = 0,1176. |iL=2,5; = c = 255"-/. 



j(SJL)=0,2375; fc(fii)=0,676. 
P(0,1176)-4-l,176-i^ + ^^^' .0,640 + 

+ ^f2T--0'135=2,93; ^=0,52. 

ror s = ~ = 0,069, j*=2; c=i;)5,;V7.,,. 
j(|x)=0,3069, /:{a)=0,808:>; 
;i.0,0«9)=4_«,0.J- . 0,50.+ 

2 • . /to 

for « =..^ =0,0377; (*=1,75; '•,=V c = 144,67.,,; 

y(a)=0,3578; /c(a)=0,8925; 
P(0,0377)=4— 0,:577-fO,010==3,6:5; 
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It is thus seen that the coefficient of the formula 
for the pressure drops sharply with decrease in velocity 
near the li^iiting velocities; after which it drops more 
slowly. Thus, when the velocity decreases from 333 me- 
ters per second to 318 meters per second, the coefficient 
drops frora 0.58 to 0.56 - decreasing hy 0.02 that is, 
3*4 percent of its value. The same numerical value for 
the drop is ohtained on changing the velocity from 196 
meters per second to 145 metors per second', although the 
difference of these velocities is 3,4 times greater. It 
may be noted, moreover, that at a velocity of 145 meters 
per second the pressure coefficie-nt is already near the 
value v/hich is computed "by the formula of Kirchhoff 

— - — =: 0.44; the differences of these values is equal to 

77+4 

0,024, abo\it 5 percent of the greater of them. - The total 
increae:it of the coefficient 0.58 - 0.44 = 0.14 is ahout 
33 percent of its lov/er value. 

Thus , at not i^_oo large ve^locit ies the coefficient in 

sur e on a '^l^te , or what is eauiv- 
9:i£2:i J Sl^sis. t a_nc^ of a ica £e ous J^e d ium t__o t^he mo t^i on in 

k21i^.liiJ^2:S. the Jiesisiance of the medium follows, appr oxi— 

ih^ s(i!i§:^.^ Isii- iih§,^ ^ii© veloci.t;y; of motion _of the 
is near the velocity of s ound , h^v/e ver , t he res is t- 
A?icrecise_s _in a ver^; parked iijanner , This conclusion 
is entirely confirmed hy the available experimental data, 
as is shown later. 

further is noted a relatively simple formula which 
for the assumed accuracy of computation gives results 
entirely agreeing with those obtained by formula (88): 



Compute P(sq) by (89) and by (88); and for comparison 

write the results one below the other. The following table 
is obtained 



P(so) = 4 ^ lOso + 7so^ 



(8S) 



So 

P(so) by (38) 
P(sJ by (89) 



0.2 



0.1818 



0.1176 



0. 069 



0.0377 



2.28 



2.42 



2.93 



3.34 



3.63 



2.28 



2.41 



2.92 



3,34 



3.63 
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The difference obtained in the value of P "by these 
fonuulas in no way affects the value of the resistance 
coefficient for the assumed approximation. 



Tlie final resistance equation now assumes the form 



E = 2loiVo^ ' (90) 



wnere 



TT + 4 - lOSg + 7Sq^ 



V 2 

^0 



where c is the vel^ocity of sound propagation at the dis- 
tant points of the f 1 ov: , equal to 333 meters per second 
•if the notion of the t)late takes :?lace in the atmosphere. 

If is not ver;- large (vq <^«^, the term 7sq^ 

is negligilDle within the limits of accuracy; then 

R . ^ 2lpaVo^ (901) 

TT + 4 - 2-^-^ . 

2he approximate formula (90) may "be obtained from the 
exact expression for the pressure of the gas jet on the 
plate in exactly the same manner as the corresponding for- 
mula (69) in the previous section. It is not 'difficult to 

show that if in the denominator of relation (87) x^ ^ is 

n , 0 

replaced approximately hy 

x^(0) + Sox'^(O) + -f! xo"(0) 



the conputation " is ■ carried out and m is set equal to 0 
formula (90) is arrived at. The difference will be only 
that the denominator will be found oqual to rr+U- 10so+ l^'^^o^ 
"but this is of .small significance for assumed accuracy of 
computation. 
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Investigation will now te made to find what formula 
(90) will yield if it is attempted to apply it in deter- 
mining air resistance to motion, using it teyond the true 
limit of its applicah ility . . 

The fundamental factor to te considered is the coef- 
ficient K 



TT + 4 - lOs.Q + 7so^ 

As long as the velocity of the plate does not exceed the 
velocit:^ of sound this coefficient increases at first very 
slowly, then much more rapidly with increase in velocity, 
as has Deen shown, the limits of its variations "being 
given hy the extreme values 0.44 and 0e58, This increase 
continues even after Sq goes heyond the value 0.2 cor- 
responding to the equation Vq = c; thus for 
V - zz ahout 500 meters per second, = 0.45, K =3 0»77; 

for s^ = ^, v^ = ^'/^c - ahout 629 meters per second; K. 

attains its maximum value 0.88, twice its value for small 
velocities. J'urther on K decreases and for Vq = 2.5c 

alDouv 833 meters per second, Sq = 1.25, there is ol)- 
tained K - 0.56. 

The aoove results qualitatively are in sufficiently 
good a::reement with test results. This is all the more 
int er GvS t ing in view of the fact that the tests were con- 
ducted under conditions very far removed from those of 
the theoretical problem considered since, in fact, K . vras 
computed from ohs ervat i ons on the flight of artillery pro- 
jectiles. The results were ohtained from hallistic tests 
Tby Zahudsky (reference 11, pp. 47-57, tabl^ 4, and fig., 30) 
and VTove mainly used for the purposes of com.parison. The 
chanj:,^e in the coefficient K for velocities not exceeding 
240 meters per second is in fact almost inappreciable; it 
then starts to grow very rapidly, increasing 2,8 times for 
a chr.nge in the velocity of the projectile from 240 to 420 
meters per second; thereafter it remains at the same level 
until the velocity exceeds 550 meters per second and then 
drops, giving for 1100 meters per second the same valiie as 
for 340 meters per second. Thus the actual change in K 
stands out with great sharpness: the law of the proportion- 
ality of the resistance to the square of the velocity 
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clearlj'- aiDplies f or. velocit ies not too large ,v/hile for 
velocities near that of sound the coefficient increases 
muchnore rapidly. This phenomenon ' might to a certain 
degree have "been foreseen. Thus in the theoretical prolD-. 
lem the tulDos of flow separate from the moving plate in ' 
only tv;o directions, while in the flight of an artillery 
shell they separate in all meridional planes* If in each 
of these planes the motion took place as in the 'case under 
consideration the resistance coefficient should vary as IC^. 
Actually the deviations are not so large. This is hecause 
the tubes of flov/ springing from the projectile spread cut 
and therefore should press less strongly on the hody of 
the projectile near its contours. Better results could 
not he expected for the reason that applied formula (90) 
is outside the limits within which its applicahility has 
"been proved. 

Acide from ■••he alDOve reasons there is yet a further 
deviating factor: namely, the viscosity of the air and 
its friction at the side of the moving tody. Ov/ing to 
"the viscosity there should he formed "behind the plate vor- 
tices v/hich lower the pressure in this region; and hence 
lead to an increase in the resistance. This condition 
already, shows up at the smaller velocities such that, as 
shown by the tests of Tihot, the coefficient. is equal 
appr orcimat ely to 0,64 instead of 0.44 if the velocity 
fluctuates within the limits of 0.5 to 11 meters per second 
(reference 11, p. 14). Jor large velocities the effect of 
the viscosity i/ould presumably be not so large. 

In concluding this part, a method is indicated for 
deriving a theoretical formula for the resistance in the 
case where the velocity of motion of the plate exceeds 
the velocity of spund. In this case, for the same reasons 
as for the case of a gas flowing out of a vessel, no con- 
tinuous steady motion should be expected. As in the 
previous case, a certain partition surface should be formed 
dividing the region of the flow into two parts in each of 
which the motion possesses a different character. This 
surf r.ce , consisting of the two sheets shov/n in figure 6, 
will bo considered as enveloping the sound waves. Within 
the compressed and heated air -layer separated by the sur- 
face from the atmosphere the motion will be steady and the 

variable t will everywhere be loss than ~- — on the 

2p + 1 

surface itself T = i and the relative velocity of 

2p + 1 
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the air particles penetrating the layer is equal to the 
sound velocity v/liich v;ould he ohscrved at t}iis place in a 
gas at rest and which is equal to the velocity of the mov— . 
ing plate. On passing through the "boundary of the layer 
into the outer at rr.ospher e , a sharp drop in pressure .is 
encountered; here the motion will he unsteady. Under the 
same assuraption as in the case of the outflcv/ from a vessel 
it is found that the angle at which the flow tuhes in their 
relative uotion (for stationary plate) intersect the ocund- 
ary of the region of the condensed layer v;ill he constant 
at all points of the partition surface. This additional 
condition is sufficient for a mathematical analysis of the • 
motion within the separating air layer and, therefore, 
also for the solution of the prohlem of the air pressure 
on the plate. It may he remarked that the very existence 
of the partition surface and condensed air layer are "by 
no means to he considered as hypothetical, since the ex- 
istence of these phenomena has heen firmljr estahlished hy 
Mach and other careful investigators. 
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PART V 



APPHOXIKATION kETEOD OF SOLUTION OF GrAS J3T PSOBLEKS 



If the velocities of the gas flovr are sufficiently 
below the limiting velocity c defernined "by the equation 

T = — \ — the solution formulas may be presented, approx— 
^ 2p+l 

imately, in a more simple and compact form by introducing 
a certain complex variable. 

In part I the folloviing fundamental equations connect- 
ing the derivatives of the velocity potential and stream 
function with respect to the independent variables T 
and 6 were derived: 

^ = 2T(i-.Tr^ ^ 

i3T .2T(l-T) S9 



These are formulas (11) of part I. The follow- 
ing notation is introduced: 




(91) 



where '^^ is the maximum value of corresponding to 

the boundary of the jet^ The preceding formulas then can 
be expressed by 

be ha 

bcp ^ b^V 
— ^ c K — 

ba b9 
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where 



K = 



l^(2p+l)T 




In the region of flovr the coefficient. K varies; 
"but j,f jbhe veloc j-t i_e5. of inojtj,.gn of jbhe ^as are nojb near 
Ih^ 1 imit ing velocity K j._s conf ined wljthi.n ZILZ narrow 

rirst it is shown that K decreases vrith increasing 
S'or this purpose the derivative dK/dT is obtained: 



It is clear that the minus sign is retained, what- 
ever the positive value T, so that the foregoing state^ 
ment is correct, IJext the values of K for the extreme 
values of t admissitle in the problem considered is 
computed. Then on the hasis of previously mentioned data 
E vjill "be included "betx^een the "boundary values thus 
obtsined. 

It is necessary to proceed, for convenience of the 
computation, from the variaole t to the varialjle 
s ?= T/d—T*) to obtain, for K, the value 



dT (l^T)^P-^^ 



K = (l--23s) (1+s) 



or if P is 



set, as "before, equal to 2,5 



E = (l-5s) (1+s) 



5 



whence is obtaihed 
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So = — = 1 > K > 0.990 



1 


1 


16 P 


40 


1 


1 


18P 


45' 


1 


1 


24P 


60 


1 


_1_ 






32p 


80 




_1_ 


5 op 


125 


1 


. 1 



1 > I > 0.9920 



« := ~ 1 > K > 0,9957 



1 > K > 0,9976 



1 > K > 0,9992 



s « « 1 > K > 0,99995 

^ 200P 500 



The corresponding values of the maximum velocity v^ 
are determined by the formula v^® 5c^s^^, where, c* is 
the velocity of sound for the physical state of the gas 
at the boundaries of the jet, For the preceding values 
of Soi if it is assumed that near the boundaries of the 
jet mean atmospheric conditions prevail, the following 
values are obtained for v^ ; 



135; 126; 118; 111; 96; 83; 66,6; 33,3 meters per 

second (the figures are rounded for simplicity; c is 
assumed eq.ual to 333 m/sec). 

The approximation which is now made consists in taking 
K equal to unity, ■ Then there is obtained: 

ba- ' he 

S0 oa 
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and, therefore 



+ ^ 3P(crt iO) (92) 



where a is determined by formula (91)^ 

It is well to consider again the corresponding motion 
of an incompressible liquid for the same boundary condi- 
tions, together with the problem on the gas flow; plane 
bounding walls , flow extends infinitely in certain di- 
rections, behind the walls at which the flowing mass 
separates the pressure is constant and the liquid or the 
gas is at rest. The problem for the case of the incom-^ 
pressible liquid is solved by the relations: 



w^ « cpi + i^^l « + ifl) . 6 + ie « Ig 

dw^ 

where d « Ig 8 is the same anglp of the velocity 

V 

♦ 

with the X axis as in the gas-flow problem^ Over all 
the boundaries has some constant value; at the 

at the bounding walls 9 constant; and at the jet sur-* 
face the velocity ' v =: constant - v^; and therefore 
- 0» These are the conditions imposed on w^^ , a 

function of the complex variable ^ + 16* The method of 
obtaining such function is given by Joukowsfcy in refer--, 
ence 2, 

It is clear tha": eJ'te:: the function P is found, 
which solves the givsu problem on the incompressible 
liquid, the required solution of the same problem on the 
gas motion is obtained by setting 



cp + ii|/ « l(a + i9) 

that is, simply replacing ^ by a; then when the variables 
6 and o pass around the boundaries of their region, ^1/ 
will receive the same constant values as "^^i where d = o 
of course a = 0, ^and, therefore, T » T^^ 

After the function cp + i\|/ is found as a function of 
a + 16 the coordinates easily can be found as functions of 
the variables a and 9, the contours of the jet investi- 
gated can be obtained and the constants characteristic of 
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the problem determln«d: namely, the quantity of ^as carried 
"by the Jet and the resultant force on the platee. 

To determine', the dependence of the coordinates on a 
and 9 T is expressed in terms of a, 7rom formula (91) 

&nd since It is as.s'uiiied that. 

A=|l-(.i.3+l)T] :(l-T)»^+'=l, 

there is obtained 

Setting T » T- yields 

~ 2~— — yi^-* ^^^^ 

Turn now to formulas (?•) and (s), part I, that gives 
the derlyatives of the coordinates with respect to 9 and ^; 

djo cosO dy siiiO 

iix _ siuO( 1 —i)-^ dy _ cesO(l--:r^ 



•« 



whence, by taking into account (93) and setting x + ly ■ 2, 
9 + i* * w, 9 - i* » w« there is obtained: 

...... 2/25| = r,^^:^-4-^-'-^i^; ; • 

She iiitegration now may be carried out, since 

»r==/(c-J- <0). <r'==/", (c- ;fJ), 

where f^ is the function conjugate to f . 
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Consider, for example, the approximate solution of the 
prolDleiTi of the pressure exerted on a plate by an infinite 
gas flov/, or otherwise expressed, of the resistance of a 
gas medium to the motion of a plate. It is assumed here 
that the direction of the flow forms a certain angle X 
with the normal to -the plate. By making use of the Joulcov/— 
sky method the following solution of this problem for the 
incompressible liq^uid is readily obtained; . 

= sin\ + sin(e^id) 

u 



The regions of flow correspond, in this case, to the 
upper half plane of the variable u. The boundaries * = 0, 
which determine the streamlines CA and • BD, correspond 
to the segments of the real axis of the u region from 
u = + ^ to us=l — sin\ and from u = - 1 — sin\ to 
u = --00 the point u ± oo gives 5=0, e = ^ \. The 
part of the plate where 9 = tt/2 corresponds to the segment 
of the real axis of the u plane bounded by the points 
u = 0 and u =: 1 ^ sin\; finally, at the boundary OB 
Q = — Tr/3 and u varies from 0 to —1 — sin\^ 

On the other hand, for u real = 0 and cp^ 

varies from 0 to + co as u varies from 0 to ±co^ The 
imaginary axis of the u region likevfise gives 0; 
cp^ increases from — ^ to 0, while u runs through the 

values from + i to 0; 5 and 9 vary correspondingly 
within the limits 0 and + c^^ ^ \ and 0, It is clear, 
from this, that the imaginary axis of the u half plane 
corresponds in the plane of flow to the streamline, EO 
branching at the plate into OAC and OBD, 

On the basis of the foregoing rule, the solution of 
this problem of the gas flow is obtained by setting 

cp + i>i/ = ku^, « 5in\ + sin(9-^ia) (95) 

The expressions for the coordinates in terms of a 
and 9 now will be sought. Turning for this purpose to 
the last of formulas (94) 
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Setting, for 'briefness, under the Integral signs 

and substituting the expressions for w and w« from (95) 
gives 

- 4/. C7,cos A J (^j--7q--7r7' ■ 
Integration yieldr 

. . o- *<'cos2A-|-sinA . ^ <4-isinA 

J <'4-vsiirA , , 

— 2tkcosAC.,arcf(f — S — ^ hi. 

Or by substitution of the sum and difference of the 
arci-aiigents in the foregoing equation and multiplication of 
tie entire equation by i: 

k [siiiA+sin(6-w)J-^^-''* [8inA4 8in(0Ti^ 

^ cos2>.-isin/^-'-'^ ... cos2' a— tsinV-'"'^ . 

- 8inX-hsiD(0-/ ;r ~ ^^'"'"'^ sinX+8in(64-w, + 

• (90) 

./^ . ^ f (e'-fe-V*+2;8inA cos/. J , 
^ V n » lco82A-fe-"-|-»8mA(c'+e ')«'• » I 

It It not difficult to thov that both arctangents enter- 
ing the foregoing equation ereryvhere rary contlnuouelyt that 
as a approaches »^ vhaterer the ralue of 8, She first 
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TT TT 

of them approaches "g ^ second g ^» ^'^^ 

therefore, L (the arbi trary constant of iategrat ion) 
may be determined so that for a=co, '2=0. 

It is noted further that the first of these arctan- 
gents iiov/here attains the value tt/2, since the denom-- 
inator of its argument nowhere "becomes zero; the second 
passes through tt/2 on the curve defined "by the eq.uation 

Ssine + sin\(e^+e""^) = 0 



!?he expression for the length of the plate 21 is 
now set up. For this purpose, by formula (96) there is 
determined 2l = ^Iz^-i- iz^ equal to the difference in the 
results of the substitution in the expression — iz of the 
valuer, a = 0, 9 = tt/2 and a = 0, 9 - - Tr/2. It may be 
noted that on the basis of v^rhat has been said of the vari- 
ation of the second arctangent of formula (P6) 



I 




( e — e ) e cos\ 

=; XT 



i~ie +sinX(e +e )e 



=.=0. 



The same substitution in the first of the arctangents 
gives zero as a result, Wtth this in mind, it is found, 
after simple reduction that: 

21 72a C, + - 
V = 4^ _i d + TTCosX ^ 

k 2 2 

whence, from formulas (93>) 



21 Jj^llo « 4 + .Tf.ccs\(l-T^) ^ 



Turn now to the computation of the resultant force E 
on the plate^ For this purpose use is made of the formula 
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P = Po(l~T) 

where is the pressure at the critical point; it is 

determined in terms of p^^ the pressure behind the plate 
and prevailing over the entire gas medium at rest "by the 
formula 

p^d-T^) = Px • 

Por deteriaining E 
1, 0 



E 



= J p^Cl-T)^"^* dy +J Po(l-T)^'*"^dy . . -2pil 

94 • e=-5 



where \^ and \^ denote the corresponding lengths of. 
the parts of the plate OA and OB from the critical point 
to the ends. 

By carrying out the integration hy parts there is 
obtained 

H = Vjl-r^'^ (la + lJ-Sp.l 
> P 

+ (P+l)Po / dT(y - y ) 

Suhstitute, in this expression, the variable a; from 
the relation between T and a (foraiulas (91) and (93)), 



3 « 
(1-t) dT « - 2Tda = dcr 



(C^^e + C.e^^) 



and the limits of integration v;ith respect to a are 

and 0; moreover it may be noted that the first two terms 
in the fgrmula for R cancel, since + = 21, 
Po^^-'^o) = Pi* Hence 
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or, integrating by parts, 

CM 



^M _ // -,i s 1_ I 



Substitution gives the result 2l/(Ci+Cg). 

With regard to the remaining integral, by equatioa (94), 
If 6 a tt/2 and It is remembered here that * = 0 



For 8 « — tt/2, >|/ = 0 there Is obtained 



Making use of these formulas and integrating again 
by parts, in the expression for R, reduces it to the form 



2 2> 



oo 



2' 

From equations (95) there is 
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and, therefore, by substitution, and introduction in the 
integral of the variable | = '"^t there ia obtained: 

1 

0 

These definite Integrals haye, respedt ively, the values 

1 f sinX— cos2 X . 1 /r. X\ ) 
2c6^A 2-f2sinX cosAVi 2)] 



and 



1 I sinX-f-cos2A • - 1 J_/"_l^M 
2ws^ i 2-T-2sinX + cosX \i^2)i ' 

and' hence, on adding, give 

2C0S-X \ ~ 2cosX/ ' 
SubstitutlQn of this expression in the formula for R yields 



Since 

2 
"0 



C,H-C2=-==, 2/.V^=4A;+riM»sX(l— tJ-^ (cm. (97)), 

l+(l->o) -^ 



therefore 

7/ = ^ Acos/V^(P+l)2?.,, (98) 
7rom the formula Just glTen for 2\i 
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Moreover 

Since from the definition of the constant a (See pt, I.) 

^ ^ KYp 0 

V-1 



and the density at the jet surface , equal to the 

density of the distant regions of the flov;, is connected 

with the density at the critical point by the foruula 

Pi " Po^^^'^o) 



Finally "by taking into account the equation ^ SaT^ == v^, 

the velocity at the jet surface and at the infinitely dis- 
tant -ooints of the moving gas mass, there is obtained from 
(98): 

TTC0S\ 

E ^ SlVo^Pi 

4(1-Tq) 4-ttcosX 



This formula for p = 0 passes over, as it should, 
into the formula of Lord Sayleigh for the flow of an in— 

compressi'ble liquid, and for x 0 gives the approximate 
solution of the problem of the pressure of a symmetrical gas 

TT C O S X. 

flow on a plate. Computing the coefficient ^ — , — 

■ 4(l--T^)^+tTC0S\ 

.for the values of Sq assumed at the "beginning of this 
part and for which this approximate method is applicable 
and considering only the case of symmetrical flovr leads to 
the following result: lor a change in s^ from 0 to 

"■gp and flow velocity from 0 to 136 meters per second 

the coefficient ~ — ^ — — fluctuates within the limits 

4(1-T^)^+Tr . ■ 



108 



NACA ^TM No, 106:^ 



12. 9,tAIL9,i increase at the end "being greater 

than at the beginning. Thus for = -~ and v = 96 

° 243. 

meters per second it still equals only 0.449, Hence for 
not very large velocities the law of the proportionality 
of the resistance to the square of the velocity is found 
to he almost exact* 



SUPPLEKSNTAHy SEKAHKS 
1. 



Part II: It is of interest to note that the function 
y will always have real roots within the limits of the 

1 

variation of from its critical value to 1, pro-. 

2p+l 

vidod n is sufficiently large. Thus for functions with 
integral n it is true for n > 1, The numher of roots 
increases infinitely vrith n. These results are ohtained 
from Pcrtor^s article (reference 12), It is readily con- 
cluded that the solution of the pro'blems on the flow of a 
gas out of a vessel and the resistance of a moving pl?.te 
in air, given in parts III and IV, are not applicable out- 
side of the limits indicated in this paper "because of the 
divergence of the series e:^pressing the stream function 
and velocity potential. 



2. 

Part V: The expression K =[l - ( 28+1)1] ( I^t)**^^"^ 
which, in presenting the "approximate method," was. accepted 
-as equal to unity actually will be eqtial to unity in two 
cases: 

1, If 3 5= 0; This is the case of the motion of an 
incompressible liquid, since the formula for the density 
P ^ Po^^'^'^)^ reduces to the equation p = constant, 

2. If 3 = - In this case p = kp ^ ^ k/p. If 

the moving matter is an ideal gas, then In order that this 
condition may be satisfied, it is necessary, in some manner 
to re;iiove the heat from the flowing mass of gas. To create 
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such a state of motion of the gas is not actually possible. 

The prohlem, however, arrived at on setting p = — ~ is of 
interest from q^uite another viewpoint^ 

The initial equations (7) of part I is considered. By 



substitution in them for t its value 




2a equal to —1 is taken and, for "briefness, the derivatives 
of cp are denoted v/ith respect to x and y by p and q 
respectively. The equations then become 



P ^ c^i/ ^ ^ ^ 

yl+p^^p ^ yi^pTq^ 

pdy— qdx 
y'l+p^+q^ 

.H.ence it is clear that if we put cp = z , then 
y, 2 v/ill he the rectangular coordinates of the points 
of a minimal surface. 

Formulas (91) to (95) of part V, on substituting -u 
for lead to the following relations: 




Khere the arbitrary constants are given soiiaewhat different 
valu.es from those in the formulas of 

If a+ i9 = t, a— i9"t^, then 

2 + i\|; = f(t) 
x-f iy=: / e'^'^ f J(t)dt + Cg / e^'^^f ^ t ( t^ )d t ^ 
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where 

q . 

u « + q^, 9 = arctg-- 

For the square of the linear element of the surface 
there is found the expression 



If the xjr plane is horizontal the curves z = constant 
will l)e the horizontals of the surface; ^^ = constant are 
their orthogonal trajectories. 

From the foregoing equations minimal surfaces of various 
shapes may be derived. 

1, Setting 

f(t) = 

yields 

n+i • llt=l 

X + iy = f£i (z+iAl/) 4. ilEs (z^i\l/) 
n+1 n— 1 

. For n rational various shapes of algebraic surfaces 
are thus obtained. An exception is the case n = 1, the 
surface then "being transcendental. 

Setting f(t) = At., gives for- real A the catenoid 
and forAthehelicoid. 

S, A second group of minimal surfaces obtained from 
the above formulas is of much greater interest, With the 
aid of the latter the minimal surface described within a 
certain given polygonal contour may be sought. The latter 
should consist of horizontal and vertical straight segments 
(the xy plane as before is taken to be. some horizontal 
plane). On setting for simplicity Uq = and hence 
^1 — Cg ?5 1/2 in the above formulas the following is 
noted: On each horizontal segment of the boundary contour 
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there will' evidently be. 2 = constant and 9= constant; 
:.if, hbv/ever,' the segment under consideration is vertical 
then on it ^l/ = constant, + = ^- and a = 0^ Thus 

the regions z + i\l/ and a + i9 will te bounded by 
straight lines. By finding the conformal transformation 
of these regions on the upper half plane of the auxiliary 
complex variable by the known method, the problem 

to the effecting of quadratures is reduced. As a very" 
simple example, the surface described in a pentagonal 
contour of the following 5hape l.s obtained: one of its 
sides is the segment of the y axis bisected by the origin 
of coordinates; from the end of this segment are drawn 
two equal sides parallel to the z axis; from the ends of 
the latter two infinite lines parallel to the x axis are 
drawn thus completing the contour^ This surface is ex- 
pressed by the follov/ing equations; 



s in 



/ ^ o -i^ J ^2 



sinh 



Sxv'l^k^ ~ ^ a " a ' a 



a a 



. i z , i z ^l/ 
2kisn — -dn — cn ^ 



1 - k^sn^ — sn — 
a a 



They are readily obtained with the aid of the pre- 
ceding general formulas if 

a + iG = 

y 1-k^ 



2 + i 



ai 



ds 



^/ (l--s")7l^k^s^) 



In conclusion, it may be noted that the given con—* 
ditions for the Surface may be somewhat varied. Thus, 
among the conditions, the requirement, that one of the 
horizontals be a line of curvature of the surface may be 
included,. The plane of this horizontal will then intersect 
the required -surface at a constant angle and, therefore, it 
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vlil te found that, for a certain given value of 2, 

+ q = constant and hence a « constant. In exactly 
the same vmy, if it is known that one of the curves,- 
11/ = constant, is a plane curve, then, as is easily shown, 
along this curve the angle 9 will be constant. 

Translation by S. Eeiss, • 
National Advisory Committee \ '.' 

for Aeronaut ics. 
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Figure 2. 
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Figs. 6,7 




